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A METHOD FOR THE CALCULATION OF LATTICE ENERGIES OF COMPLEX 
CRYSTALS WITH APPLICATION TO THE OXIDES OF MOLYBDENUM 

By William S. C hanky 


SUMMARY 

A theoretical study has bee a made o f molybdenum 
dioxide and molybde num trios id e in order to extend 
the knowledge of conditions and factors invoiced in 
the oxidation of molybdenum . The lattice energies 

hare been calculated by using an improved and more 
generalized method which is based on the electro- 
static model. The crystal structure has been exam- 
ined and the details of the structure have been cor- 
related with the valence mechanisms contributing to 
the lattice energy. 

There are several energy-producing valence mecha- 
nisms which contribute to the lattice energy. / hose 
studied herein are coufombic energy , polarization 
energy , Van der Wauls energy , and repulsion energy. 
The methods used to compute these energies have 
been modified and improved making them applicable 
to many types of crystals which could not be readily 
studied, with the old methods. In particular, the 
new method permits the calculation of lattice energies 
of distorted crystal structures such as those found in 
the ox ides o f m olybde n um . 

The data reported by X-ray crysta/lographers have 
been transformed into a vector system. This system 
is easier to apply in the energy calculations and 
easier to use in sketching the crystal structure . 

('ova lent effects have been examined in the molyb- 
denum oxide crystals but a quantitative calculation 
of covalent energy was not attempted. A postulated 
covalently induced electrostatic energy was briefly 
investigated and appears to be feasible. 

It is felt that the overall study has provided much 
insight as to the nature of the bonding mechanisms 
responsible for the stability of these oxides. In 
particular , the layer structure and weak interlayer 
bond found in molybdenum trioxide are believed to 


be responsible for the high vapor pressure of tins 
oxide. 

Perhaps of mote general significance , the methods 
developed in this research should be useful in extend- 
ing the knowledge of complex crystals which have 
not been subject to this type of study be fore. 

INTRODUCTION 

The electrostatic theories of valence have been 
used for* many years to compute lattice energies 
and to study the relationships among crystal 
structure, lattice energies, and physical properties 
of crystals. The established methods are limited, 
however, to relatively simple crystal structures 
possessing great symmetry. The mathematical 
methods that are employed take advantage of the 
crystal symmetry to express the energy relation- 
ships in terms of mathematical series which may 
he resolved by established techniques. Although 
the complexity of the established methods limits 
their application to relatively simple symmetrical 
crystal structures, many real crystals are amenable 
to these methods. 

In this investigation, the electrostatic theories 
of valence were applied to the dioxide and trioxide 
of molybdenum. Because of the unusual and 
distorted crystal structure of these oxides, the 
usual method of computing the components of 
lattice energy could not be applied. A new 
method has been developed which is more general 
in its scope. This method is discussed in great 
detail and it is applied to these two oxides of 
molybdenum. Because of the more general char- 
acter of this method, it should he possible to 
apply the electrostatic theories of valence to the 
many complex crystals which have not been 
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investigated by this means before. 

Molybdenum is a metal which lias attracted 
interest. It possesses properties which make it 
a potentially useful structural material in high- 
temperature environments. However, in an 
environment with very high temperatures and an 
oxygen containing atmosphere, such as the envi- 
ronment for very high-speed terrestrial flight, 
molybdenum combines with oxygen catastro- 
phically. 

The oxidation of molybdenum has been studied 
by many investigator's using empirical and experi- 
mental methods. Apparently no theoretical study 
has been made of this problem. A theoretical 
studv was undertaken here to extend the knowl- 
edge of the conditions and factors involved in the 
oxidation of molybdenum. 

The author wishes to express his appreciation to 
Dr. (). K. Rice, Professor of Chemistry, Univer- 
sity of North Carolina, for constructive criticism 
and helpful advice. 

SYMBOLS 

A m M adehing constant 

a 0 Bohr radius 

a,b,e crystal cell dimensions 
B repulsion constant 

c conversion factor, keal/erg 

I) vector between two ions 

d distance between ions 

bh ionization potential 

e electronic charge 

F electric field vector 

Ci position vector of a polarizable ion 

g repulsion exponential law parameter 

— > 

T general position vector of an ion 

^,j,k unit vectors in an orthogonal system 
J Van del* Waals dipole-dipole coefficient 

Ij Van der Waals dipole-quadrupole coeffi- 

cient 

M u number of bonds between ions i and j 

Mo' molybdenum ion position vector 

X number of polarizable elections 

Nat Avogudro's number or number of polar- 

izable ions 

n repulsion power law parameter 

O’t oxygen ion position vector 

P number of outer shell electrons in an ion 


p r repulsion coefficient defined by equation 

( 74 ) 

q electric charge 

R interionic distance at equilibrium 

R vector from the origin of one cell to the 

origin of another cell 
Rj ionic refraction of the jth ion 

r ionic radius 

V energy 

l T c coulombie energy 

r <M Van der Waals dipole-dipole energy 
C <I(J Van der Waals dipole-quadrupole energy 

r„ polarizat ion energy 

C r repulsion energy 

C s stabilization energy 

V t total lattice energy 

x,y,z Cartesian coordinate system axes 

Z valence 

a polarizability 

«,d, 7 mteraxial angles of the crystal cell 

a, ( 3 . 7 unit vectors in the directions of the 

crystal-cell axes 

ri parameter, layer number 

p scaling parameter 

t number of ions per unit cell 

p electric dipole moment 

£ number of moles or stoichiometric units 

per unit cell 

THE OXIDATION PROBLEM 

In general, all metals except certain noble metals 
have i natural tendency to react with oxygen gas 
and form an oxide. The usefulness of a metal as a 
structural material is greatly dependent, on the 
rate it which this oxidation process proceeds. 
The basic method of controlling the oxidation proc- 
ess is o interpose a harrier between the metal sur- 
face imd the atmosphere. This barrier may be an 
nrtilicial device, such as a coat of paint ; however, 
for iretals known to he inherently oxidation re- 
sistant, the barrier is usually a tough tightly 
adhei cut oxide film. 

Th 4 oxide of a metal is a very important factor 
in tin oxidation characteristics of the metal. The 
oxide film may be loose, porous, or flaky and pro- 
vide little protection to the underlying metal. A 
protective oxide will form a tough and tightly 
adherent film. If a break in this film should 
occur then subsequent oxidation will repair tin* 
defee . This is said to be a self-regenerative film. 
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The oxidation process that occurs in the presence 
of a self-regenerative protective oxide film is 
controlled by the rates of diffusion of either the 
oxygen atoms or the metal atoms or both. The 
force that drives the diffusion of atoms through 
the film is related to the gradient of the chemical 
potential of these atoms across the film. As the 
film becomes thick and the gradient of t lie chem- 
ical potential small, the oxidation process may 
reach a virtual standstill. A on protective oxides 
may occur where the oxide is in the gas state 
rather than the solid state. Molybdenum tri- 
oxide exists as a solid at lower temperatures but 
readily vaporizes at higher temperatures affording 
practically no protection for the metal. At high 
temperatures, molybdenum oxidizes very rapidly 
and this imposes a serious limit on the usefulness 
of the metal as a structural material. 

Tli ere are several oxides of molybdenum known 
to form under various conditions with compositions 
ranging from a metallic oxide Mo 3 () to the tri- 
oxide Mo0 3 . The formation of oxides on the 
metal surface is probably a complex transition 
through several of the possible compositions. At 
high temperatures the dioxide changes to the tri- 
oxide. The trioxide is unstable in the sense that 
it does not remain solid and certainly does not 
form a protective coating. The trioxidc melts 
at a temperature of 1,0()N° K under pressures as 
low as 10 mm Hg. At a temperature of 1 ,42N° K 
and a pressure of 1 atmosphere, the molten molyb- 
denum trioxide vaporizes. Direct sublimation 
occurs at 0.3 mm Hg and 973° K. 

There are two general areas of importance in 
research on oxides as they are related to the 
oxidation of metals. The first pertains to the 
physical nature of the oxide itself and to the 
underlying reasons for the failure of the oxide to 
form a protective film. The second pertains to 
modifications of the oxide to make it more stable 
and to reduce the diffusion controlled oxidation 
rate. The feasibility of research in this second 
area has been demonstrated by Wagner in his 
theory of oxidation (discussed in ref. 1) where it 
is shown that small additions to the base metal 
can modify a diffusion controlled oxidation rate. 
Also, some double oxides are more stable than the 
single oxides forming them. These two areas of re- 
search are related in that the first provides a 
foundation of knowledge on which the second 
may lx* bast'd . 


The present research is devoted to an explora- 
tion of the kinds of forces holding the oxides of 
molybdenum together, of the relative magnitudes 
of the energies involved, and of the relationship 
of these features to the crystal structure. 

INTRODUCTORY DISCUSSION OF THEORETICAL 
CONCEPTS 

An introductory discussion of the underlying 
concepts which form the basis for the theoretical 
developments in this investigation is given here. 
Readers who are familiar with this field of work 
may desire to omit reading this section. 

The electrostatic theory of valence is a very 
old theory and contains contributions from many 
early workers. The method of computing the 
bonding energy of molecules and the lattice en- 
ergy of crystals has long been established for cer- 
tain compounds. This method is capable of 
yielding accurate quantitative results in many 
cases, especially for crystals that, are essentially 
ionic in the nature of t heir bonds. Jt is possible 
to account for the structural configuration of 
many presumably covalent compounds on the 
basis of the electrostatic theory. 

The model used in the electrostatic theory is 
based on the assumption that ions are formed by 
the transfer of electrons from the more electro- 
positive element to the more electronegative ele- 
ment, thus forming cations and anions, respectively. 
The ions are assumed to exist in the crystal 
as quasi-spheres with a moderate amount of 
hardness and with a variation in size over a small 
range. Except when the ions are polarized, the 
electric charge is assumed to be concentrated at 
the center of the sphere. In the presence of an 
electric field, the ions may he polarized by ac- 
quiring an induced separation of the centers of 
positive and negative charge. Tin* ions contain 
“orbiting” electrons which, on a short-time basis, 
act with the nucleus to produce a revolving elec- 
tric dipole and in two adjacent ions the revolution 
of these dipoles may become synchronized to 
produce an attractive force between the ions. 
When forced together the ions repel one another 
by repulsion forces (other than coulomb forces) 
which presumably act in compliance with the 
Pauli exclusion principle. In most computations 
the temperature is presumed to lx* 0° K. The 
ions in crystals vibrate about their mean positions 
with zero-point- vibrational energy. In cases 
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where t lie cation has partially filled d-orbituls, 
crystal field effects may he present. The crystal 
field effects are not present in the trioxide of 
molybdenum and are assumed to he ineffective in 
the case of the dioxide. 

in lattice energy calculations based on the elec- 
trostatic theory, there are several energy contri- 
butions which are computed individually. These 
include eoulombic energy, polarization energy, 
Van der Waals energy, and repulsion energy. 
Zero-point vibrational energy and crystal field ef- 
fects are not included in this work. The theoreti- 
cal concepts are developed in the text and the 
application of these concepts is carried out in 
appendixes A to D. 

COULOMBIC ENERGY 

The eoulombic energy is the electrostatic po- 
tential energy of a molecule or stoichiometric unit 
located in an infinitely large crystal. Since the 
charge of each ion is assumed to be concentrated 
at the center point of that ion and potential energy 
contributions of neighbor's are additive, the com- 
putation of eoulombic energy amounts to a sum- 
mation of potential energy contributions from an 
infinite array of point charges. The electrostatic 
potential energy of one ion due to another is 
proportional to the reciprocal of the distance be- 
tween them and decreases in magnitude as the 
distance increases. However, for one reference 
ion in an infinite crystal, the number of neighbors 
increases rapidly with the distance away from that 
ion. The overall crystal is electrically neutral and 
there are positive and negative energy contribu- 
tions resulting from t lie different pairings of pos- 
itive and negative ions. The net effect is that for 
a reference ion the electrostatic potential energy 
increases in magnitude as neighbors are added in 
an array about it and rapidly levels off to a fixed 
value. In simple crystals, the potential energy 
contributions may he arranged in the form of a 
mathematical infinite series and the limiting 
value of the sum determined by analytical tech- 
niques. This method lias been used, in several 
variations, to effect such a summation for most 
common types of simple crystals. 

Where a crystal is found to be distorted with 
very little periodicity of the ion spacing within a 
single cell, the labor involved in the method just- 
described would be very great and the results 
would he applicable to only one anomalous 


crystal. In this research a method of effecting 
the summation is developed that requires no 
symmetry at all within a crystal cell. This method 
takes- advantage of the periodic arrangement of 
the colls themselves and the electrostatic poten- 
tial energy of a cell of ions in an infinite crystal of 
cells is determined. The summation of potential 
energy contributions for a single reference cell as 
neighboring cells are added in layers around it is 
expressed in terms of a distance parameter rj. 
The calculation is made on an automatic computer 
with the lattice expanding until eovergenee is 
established. The final contribution of long- 
distance neighbors is determined by extrapolating 
to the limit of the potential energy as 77 approaches 
infinity. This method is applicable to anv 
crystal and is limited only by t lie computational 
speed of the computer. 

VECTOR SYSTEM OF EXPRESSING CRYSTAL STRUCTURE 

( Y nsidcration of the problem of computing the 
eoulombic energy of a crystal should make it 
apparent that a system is needed to express the 
position of each ion within a unit cell and to ex- 
press the position of any cell with respect to a selec- 
ted reference cell. A vector system has been selec- 
ted because of the convenience of the vector notation 
and established vector methods of computation. 
The system adopted is actually a compound 
system consisting of reference systems located at 
each {joint in a larger reference system. Each 
cell in the crystal constitutes a reference system 
with the cell edges parallel to the cell reference sys- 
tem ; xes. The origin of the cell is the origin of the 
cell reference system. The position of each ion 
in the cell is specified by a position vector in the 
cell reference system. In the larger reference 
system, a single cell is selected as a reference cell 
and 1 he origin of the cell is the origin of the svstem. 
The axes are again parallel to the cell axes but 
distances along each axis art* measured in units 
of e dl dimensions. In this system, a vector 
specifies the position of the origin of a cell with 
respect to the origin of the reference cell. 

T 1 e advantage in using this compound system 
is tint the relative position of any one ion with 
respict to any other ion may be simply expressed 
as an algebraic sum of vectors. The notation is 
particularly convenient for the methods of corn- 
put a don developed in this research. 
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The position of each ion or atom in the cell of 
a crystal has been determined, for many sub- 
stances, by experimental measurements using 
X-ray diffraction techniques. The experimental 
data which are used in this work are given in 
appendix A. 

POLARIZATION ENERGY 

Polarization energy exists for an ion in a crystal 
when there is an electric field present at the posi- 
tion of the ion. The electric field acts to 
cause a separation of the centers of positive 
and negative charge in the crystal producing 
an electric dipole. This electric dipole then 
acquires an electrostatic potential energy in the 
electric field. The net energy effect is called the 
polarization energy. Most simple symmetrical 
crystals have a symmetric distribution of positive 
and negative charges for which the associated 
electric fields cancel at the position of each ion 
and polarization energy does not exist. The re- 
quirement of symmetry that made possible the 
calculation of eoulombie energy with mathemat- 
ical series effectively eliminated crystals in which 
polarization energy would be expected. Polari- 
zation energy is to be expected in many distorted 
crystals and isolated molecules. 

The computation of the electric field at some 
selected point in an array of point charges re- 
quires a vector summation of the electric fields 
associated with each charge since the electric* 
field is a vector quantity. This summation is 
somewhat analogous to the eoulombie energy 
problem except that convergence is much more 
rapid. 

VAN DER WAALS ENERGY 

Van der Waals energy, a kind of potential 
energy, is associated with Van der Waals forces, 
a term applied to general forces of attraction and 
repulsion of uncharged molecules. This term in- 
cludes the interaction of molecules which possess 
permanent electric dipole moments. The inter- 
action of molecules, atoms, and ions not possessing 
permanent dipole moments is the kind of Van der 
Waals energy being considered in this research. 
The associated forces are also called London 
forces. The presence of a charge on an ion does 
not prohibit the existence of Van der Waals forces; 
rather, these forces are superimposed on the cou- 
lombic forces which an* attributable to the* charges. 


A mechanical model frequently used to explain 
the Van der Waals forces between atoms is de- 
scribed by the synchronized revolving dipoles 
mentioned previously. The theory of Van der 
Waals energy is discussed in almost every text- 
book on quantum chemistry and will not be 
repeated in detail in this paper. 

Van der Wauls energy is a relatively small 
contribution to the total lattice energy of an ionic 
crystal. It can he an important factor in some 
crystals in determining which of two possible 
similar structures is most stable. 

REPULSION ENERGY 

Whenever two atoms or ions are forced together 
they tend to resist this force; that is, they tend to 
occupy a certain volume in space and resist any 
action which would reduce this volume. This 
behavior presumably is attributable to the Pauli 
exclusion principle. The methods of computing 
the energy associated with this repulsion force are 
empirical. There are two methods In general use 
known as the “inverse power law” and the 
“exponential law.” Each method employs an 
expression with two unknown parameters. One* 
parameter is eliminated by setting the derivative 
of lattice energy with respect to the interionic 
distance equal to zero at the experimentally 
determined distance. It is presumed that at Uus 
experimentally determined distance, the potential 
energy forces are all in equilibrium. The second 
unknown parameter may 1 >e determined by using 
experimentally determined compressibilities or by 
using empirical constants. The inverst* power 
law was used in this research because it was 
mathematically simpler. Both laws are valid for 
ions the potential energy forces of which are in 
equilibrium. The exponential law is considered 
by many investigators to he better where the 
forces are not in equilibrium. 

BORN-HABER CYCLE 

The validity of the model used in these calcula- 
tions of energy may be tested by comparing the 
theoretically determined lattice energy with the 
experimentally determined lattice energy. The 
Born-IIaber cycle is used to determine the 
experimental lattice energy from thermochemical 
measurements of the heats of formation. The 
heat of formation of a crystal is the energy 
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released when the crystal is formed from the 
elements in their natural state at standard condi- 
tions. The lattice energy is energy released when 
the crystal is formed from the elements in the state 
of gaseous ions. The Born -Ha her cycle allows for 
such factors as heats of vaporization and dissocia- 
tion and also for ionization potentials and electron 
affinities in computation of empirical lattice 
energy. 

unique layer structure of molybdenum trioxide 

l pon examination of the crystal structure of 
molybdenum trioxide, it becomes apparent that 
this crystal is made up of layers. In nonlayer 
crystals, the metal and nonmetal atoms are in 
alternate positions so that like atoms are not adja- 
cent to one another. Layer crystals, in this sense, 
are typified by structural layers the boundaries of 
which are defined by like atoms in the surface of 
one layer being adjacent to like atoms in the sur- 
face of the neighboring layer. Layer structures 
are believed to result from directed covalent 
bonds or large polarization of the anions. This 
same structural eoneept can be carried further to 
include crystals made up of chains or crystals of 
small individual molecular configurations wherein 
the boundaries are still defined by like adjacent 
atoms. 

The importance of layer structures in consider- 
ation of their stability characteristics lies in the 
kind of bonding between the layers. It is not 
necessary that a crystal have every bond in it 
broken in order that it lose its solid state. In the 
sense of becoming nonsolid, the crystal may fail 
by separating into small structural units. Only a 
small part of all the bonds in a crystal may break 
and, in view of this, it is worthwhile to examine 
the different kinds of bonds in a crystal to sec 
which might be the weakest. 

It is generally believed that the interlayer bond- 
ing in ionic crystals is predominantly Van der 
Wauls bonding. This conclusion is reached by 
elimination of the other types of bonds. The 
most influential type of bond that might be pres- 
ent in an ionic crystal is the coulombic bond. The 
section entitled “ Interlayer Energy” is devoted to 
an effort to make a quantitative estimate of how 
much coulombic contribution there is in the inter- 
layer bond. It is shown there that Van der Waals 
bonding is not predominant. 


covalent energy 

Gov dent energy is sometimes discussed in con- 
nection with ionic crystals as a sort of residual 
energy that represents the difference in measured 
lattice energy and lattice energy computed on the 
basis of the electrostatic model. This is not a 
good explanation of the lack of agreement in 
energy calculations although where the difference 
in energy is small the error in concept is probably 
not so objectionable as when t his difference is large. 
Generally, the eoneept of the covalent bond im- 
plies a utilization of the valence electrons that 
would make the assumptions in the electrostatic 
theory invalid. However, in the ease of molyb- 
denum dioxide, there are two d-eleetrons for the 
plus-four molybdenum ion that are not a part of 
the valence electrons so far as the electrostatic 
model is concerned. On the basis of distance 
between pairs of molybdenum ions, Pauling has 
suggested that these d-electrons interact to form 
a covalent bond (ref. 2, p. 4117). If this is true, 
then, since these d-eleetrons an* not a part of the 
valenc » electrons in tin* electrostatic model, t be- 
en ergy of this ronvalent bond would exist inde- 
pendently of the electrostatic energies. It hap- 
pens tl at tin* difference in measured and computed 
lattice energies for molybdenum dioxide is much 
great ei than a reasonable upper limit to an esti- 
mate of this covalent bond energy. Because of 
this, i is worthwhile to look further to see if 
some unaccounted energy-producing mechanism is 
presen . 

A la •ge difference in an ionic bond as compared 
with a covalent, bond lies in the electron density 
in the region between nuclei. For a pure* ionic, 
bond, tin* valence electrons are virtually com- 
pletely transferred from tin* cation to the anion 
and tl e two ions are attracted by a coulombic 
force. In this east*, on examining the electron 
density along a lint* from tin* cation nucleus to 
tilt* an on nucleus, there will be a point where the 
elect rtu density is virtually zero. In contrast, 
for a *ovalent bond the electron density is not 
zero a any point between nuclei but rather it is 
great e 1 than the electron density of free atoms 
posit icned in the same locations. This implies a 
localization of charge between nuclei rather than 
charge transfer which is compatible with the con- 
ventional concept- that in the covalent- bond there 
is a sharing of electrons. 
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lu tilt’s sod ion entitled “Covalent Effects,” it is 
assumed that the d-electrons in the molybdenum 
dioxide lattice are diffusely located between the 
nuclei of closely spaced molybdenum ions. A 
rather crude electrostatic model is then hypothe- 
sized and used to compute the coulombic energy. 
On comparing coulombic energy values computed 
with different internuclear charges with measured 
energy values, it can be determined that the 
internuclear charge which yields the correct eou- 
lombic energy is of the right order of magnitude 
to imply that this hypothesized model is feasible. 

CRYSTAL STRUCTURE 

The physical mechanisms which operate to 
produce the different contributions to the lattice 
energy and the physical properties of a crystal 
are closely related to the geometrical arrangement 
of the atoms in the crystal. It is necessary to be 
able to orient- the crystal conceptually and to 
express the position of each atom in a form which 
may be adapted to subsequent mathematical 
developments. 

It is difficult to form a mental image of a com- 
plex crystal. A device frequently used to facili- 
tate description of complex crystals is to describe 
them in terms of polyhedru. Many crystals are 
made up of octahedra or tetrahedra and the 
description of the crystal may be simplified by 
expressing relationships of, say, octahedra rather 
than atoms. The relationship between atoms is 
expressed in the description of a single octahedron. 
In both molybdenum dioxide and molybdenum 
trioxide the octahedra are formed with an oxygen 
atom at each of the six corners and the molyb- 
denum atom in the center. This arrangement, is 
shown in figure 1. The octahedron shown in this 
figure is regular. It will he seen subsequently 
that the real octahedra in these oxides are dis- 
torted. The representative octahedron is shown 
in figure 2. In schematic sketches of crystals, 
corner views and edge views of the octahedra are 
used and these are shown in figure T Molyb- 
denum trioxide is represented in figure 4 where 
figure 4(a) shows the top view of a single layer 
and figure 4(h) shows the edge view of the layers. 
Molybdenum dioxide is shown in figure 5. Of 
course, figures 4 and 5 represent only a small 
section of the crystal and the extension in all 
directions should he apparent. 



Fkji'RK 1. — Octahedron representing an arrangement of 
atoms. 


i 



Fine he 2.-— Typical schematic drawing of octahedron 
used to represent atomic configuration. 
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(a) Corner view of an octahedron. 

(b) P'dge view of an octahedron. 

Fm;i kk Hepresentative views of an octahedron. 




Ficckb 4. — Schematic structure of Mo() 3 . 


When octahedru are joined then the atoms 
located at such a joint are shared by the neigh- 
boring octahedru. Frequently the sharing of 
atoms by octahedra will define the differences iti 



Fiucre 5.-- Sciiemat ie structure of Mot > 2 . 


“kinds 7 of atoms of the same element according 
to environment. By inspection of figure 4, it may 
be seen that there are three “kinds” of oxygen 
atoms in molybdenum trioxide. The three “kinds” 
are those shared by three octahedru-, those shared 
by twc octahedra, and those that are unshared. 
The two “kinds” of oxygen atoms in molybdenum 
dioxide are not so apparent- in figure 5. In this 
case, the distortion of the octahedra is significant 
in establishing environmental “kinds” of atoms. 

The structure of crystals is established by spe- 
cialists using X-ray diffraction techniques. The 
general subject area of crystallography and X-ray 
diffraction studies is somewhat complex in nota- 
tion an 1 terminology. In order to use this infor- 
mation it is necessary to translate the description 
of the Tystal from the compact notation of the 
X-ray Tvstullogruphcr into a longer form more 
appropriate for mathematical developments. This 
form nay be called the “reconstructed lat- tire.” 
In a piaeticu-1 sense, there are several features of 
the cnstnl structure that should he expressed. 

1 he siz * and shape of the unit cell must be known. 
This it: formation is given by the type of lattice 
and the cell dimensions. Next, the position of 
each atom within the cell must be known. This 
is expressed by the space group and the atom 
position parameters. For use in subsequent 
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developments in this paper, it is convenient to 
express the position of the atoms in the coll by 
position vectors. The axes of the reference 
system for these vectors are parallel to the cell 
edges and the unit vectors are 1 angstrom unit 
long. The system is not necessarily orthogonal. 

The lattices of molybdenum dioxide and tri- 
oxide are reconstructed by using the following 
method. The outlines of the cells and the posi- 
tions of the atoms are given in figures 6 to 9. 
These sketches are not easy to interpret from a 
casual glance but rather require the coordinated 
study of different views, with attention directed 
toward atoms of particular interest. It is possible 
to locate and identify each atom and its neighbors 
and then to associate the proper position vector 
with each atom. A method is given for taking 
X-ray diffraction structural data in the form in 
which it is conventionally presented and from this 
arriving at a set of atom or ion position vectors. 
These position vectors are then used to recon- 
struct the lattice, to calculate bond lengths, and 
to calculate interbond angles. 


Structural data from X-ray diffraction studies 
include the following information : 

(1) Type of lattice 

(2) Space group 

(3) ('ell dimensions 

(4) Atom positions 

The atom positions will be given according to 
the equivalent positions that are occupied and 
the atom position parameters to be used with 
the general coordinates of the equivalent posi- 
tions. The general coordinates of the equivalent 
positions are established by the symmetry prop- 
erties of the space group. This information has 
been assembled in the International Tables for 
X-ray Crystallography (ref. 3) wherein also the 
conventional symbols of crystallography used in 
this investigation are given. The information 
necessary for the present work is given in appen- 
dix A. 

The atom position parameters may be substi- 
tuted into the general coordinates of the appro- 
priate equivalent position and the result will be 
the coordinates of the atoms in units of (*<41 
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KicniK 8.— Partial projection of MoO :i model on ab-plane. 1 inch = 2 A. 


li* the position parameters are substituted in the 
coordinates of the (4c) positions, the coordinates 
of the ions are obtained. These are 


0.232 

0.000 

0.017 

0.232 

- 0.000 

— 0.017 

0.232 

0.500 

0.483 

0.232 

0.500 

0.517 


Multipl ying the coordinates by the cell dimen- 
sions gives the position vectors as follows: 

\lo‘ = l . 2 9 5a + 0 . 000/3 T 0 . 09 53 y 
Mo" = - 1 .295^-0.0000-0.09537 
M o 1 " = - 1 . 29 5a + 2 . 42 1 Is + 2 . 7097 
Mo 1 ' — 1 ,295a + 2. 421 0+2.8997 


These coordinates are expressed in units of cell 
dimensions. They may be multiplied by the cell 
dimensions to obtain the coordinates in angstrom 
units. The cell dimensions are 

a = 5.584 A 
b = 4.842 A 
<■ = 5.608 A 


In the same' manner, the position vectors for 
the ozygen ions in Mo() 2 are found. The oxygen 
ions also occupy 4(e) positions. There are two 
“kinds” with position parameters as follows: 
Four () r in 4(e), 

x = ().l 1 y=0.21 z = 0.24 

four () n in 4(e), 


0=120.94° 


x = 0.39 


v = 0.70 


z = 0.30 
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The oxygen atom position vectors are 


Oi =().6i4a4 1.01704 1.3467 
()V = — 0.014a — 1 .0170 — 1 .3467 
()i“= — 0.6 1 4a + 3.4380 + 1.4587 
oy - o.o 14a 4- 1 .4040 4-4 . 1 50 7 
^-2.17Sa + 3.3890+l. 0827 
OlV- —2.178a — 3.3890 — 1 .0827 
Og- -2.1 78a + 5.8100+1.1227 
Oi = 2.178a-O.9680 4-4.4867 


In this notation, the position vector is labeled by 
the a omic symbol of t he element involved with a 
Roman numeral subscript designating a “kind” of 
atom and a superscript designating the particular 
equivalent position. The position vectors for 
molybdenum trioxide have been determined and 
are given in appendix A. A particular bond 
between two atoms may be expressed as a vector 
by taking the vector difference of the position 
vectors of the two atoms. The length of this bond 
is rej (lily determined by using the conventional 
dot product. 

LATTICE ENERGY 

The lattice energy is the energy released when a 
crystal is formed from elements which are in a 
state of being both gaseous and ionized. To be 
exact, this is an enthalpy since there will be 
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a change in volume hut the pressure-volume 
energy is usually very small compared with the 
internal energy and it is frequently ignored. 

The lattice energy consists of a summation of 
energy contributions from several different energy- 
producing mechanisms. In a generalized form 


this 

summation is represented by the 

following 

equation: 



l-, = lT e +Up + lTdd + Cd,+U r 

(M) 

who! 

re 


r, 

lattice energy 


r. 

coulombic energy 


u P 

polarization energy 


Cm 

dipolo-dipole Van dor Wauls energy 

^ clc) 

dipole-quad rupole Van der Waals 

energy 

Ur 

repulsion energy 



It will be shown that the repulsion energy has 
the form 

I T r = — Pr(l c + 4l r t) + t)U d( i + 8Ud(|) (4) 

and each of the other energy expressions will have 
the following forms: 


Uc = 

A M e- 

p 

c r >) 


( ^a0“ 

(«) 

Up = 

P 4 


J 

? 

(7) 

Uj„= 

L 

(«) 


In these expressions 

A m Madelung constant 

e electron charge 

C constant for one type of crystal 
a polarizability 

J dipole-dipole coefficient 

L di pole-quad ru pole coefficient 

p scaling parameter 

p r repulsion coefficient 

In the following four sections these energy 
expressions are discussed separately and equa- 
tions for computing each energy contribution to 
the lattice energy are derived. These equations 


are applied in practical calculations of the lattice 
energy of molybdenum dioxide and molybdenum 
trioxide in appendix B. In the section entitled 
“Summary of Energies” these contributions to 
the lattice energy are listed, added, and compared 
with the measured lattice energies. 

COULOMBIC ENERGY 

The coulombic energy of a crystal is the energy 
t hat is attributable to the electrostatic inter- 
action of charged ions. In the model adopted for 
these calculations, it is assumed that the atoms 
present are ionized and that their charges are 
concentrated at points at the centers of the ions. 
The locations of these center points have been 
experimentally determined and are expressed by 
the atom or* ion position vector's. Because of the 
distortion in the lattices of the molybdenum oxides, 
one of the established Madelung constants cannot 
be used and the usual methods of determining 
Madelung constants are not applicable. A general 
method of computing coulombic energy and the 
Madelung constant is developed here. This 
method requires the use of a digital computer but 
the method is simple and systematic and may be 
applied to any kind of crystal lattice. The 
coulombic energy of the dioxide and trioxide are 
determined. 

Consider a crystal made up of positive and 
negative ions. There may be ions of several 
elements present. The position of each ion within 
a unit cell is known and is given by its position 
vector. The position of one unit cell relative to a 
cell selected to represent the origin is given by a 
vector 

R = iija + n 2 /3 + n 3 7 (9) 

A ^ A * 

where ft, and y are unit vectors along the cell 
axes and n u n 2 , and n 3 are integral multiples of 
the appropriate cell dimensions. Also, 

ni=/iaA| 

n 2 = 41 ) r (10) 

n 3 — / 3 e J 

where a, b, and c are cell dimensions and l u U, 
/» = ± 1 7 ±2, ±3, .... A method of assigning 
values to i u / 2 , and / 3 is given subsequently. 

The electrostatic potential energy, in ergs, of 
one ion, say the ith ion, due to another ion, the 
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jth ion, is given by tlu* following expression: 


whore 


, _ Zi/,0* 

lj_: <h, 


(n) 


/ valance of the ion 

e electronic charge, statcoulombs 

(1 distance between the two ions, cm 

and 

<lij = ! Dij| ~ (Du • L)„) 1/2 (12) 

The term I ) u is a vector directed from the it h ion 
to the jth ion and is expressed as follows: 

Dtj-Tj+R-Ii (Id) 

> 

The vectors I, and Ij are the position vectors of 
the it h and jth ions, respectively. 

Tlic potential energy of all the ions in a given 
cell, the Oth cell, due to the ions in another cell, 
tin 1 Ktli cell, is expressed as follows (here, i 
represents the Oth (*<41 and j the Rth cell): 

>*<>*=£ u, ,=2 z f Jt ‘ 2 (, 4 ) 

ij ij ( 1 1 j 

The total potential energy of the ions in the Oth 
cell due to all the other ions present is obtained 
by summing u OH over all values of R from zero to 
infinity as 

U 0 =ll00+S **<>R (l- r >) 

n=i 

The expression for R = () is considered separately 
since all terms where R = 0 and i = j must be 
discarded to avoid dividing by zero. These would 
represent the interaction between an ion and 
itself. 

Since it is desirable for the energy to be ex- 
pressed on a pin* mole basis the final energy 
expression becomes 


The coefficient 1/2 compensates for double count- 
ing. 

In cases wliere a Madelung constant is used, 
the energy expression is given by equation (5) 
which is 




A M e“ 


p 


when 1 o is a scaling parameter. Frequently p is 
the closest interatomic distance in the crystal. 
In the present development tin* Madelung constant 
is giver as 


A 


M — 


p -Ms z ; Z) l +2 L? 1 

L ij Jr=o n = i ij <hj 

(17) 


There is no advantage in expressing a Madelung 
constant unless the crystal lattice is a frequently 
encouircred type. The molybdenum oxide lat- 
tice's ai ' unique and their Madelung constants are 
not computed. 

In expanding the coulombic energy expression, 
equation (1(>), the first summation should be made 
over all the ions i and j and then this rather large 
expression is expanded by summing over all cells. 
The summation over i and j is straight forward. A 
systematic*- method is needed to sum over R and 
extrapolate to infinity. 

Consider the expressions for djj (eqs. (12) and 
(13), respectively) 


«i, J =(iS 1J - n ,,)' 2 


D,,= i,+R-Y, 


'+ A A A 

Ji~I)i« + Jj 2 /3-r I J3 7 

(IS) 

Ii = Ina+ Ii 2 /3+Ji : i 7 

(19) 


In tin's » ('([nations the position vectors are written 

in tern s of the unit vectors of the system and 

— ► 

coeflicic nts. Tin* vector R can ho written 


n av c 0 _ x A vc* z . z.i f, ^ z,z, 

‘-'S L*T «Iii Jh-.A^iY «I|J 

(iti) 

wliere 

N’av Avojradro’s number 
£ number of stoichiometric units per cell 


R-/,aa+J 2 b0+/3<'r (20) 

- > --> - > 

Using these forms for 1 ,, 1,, and R gives 

I)ij = (I it + /iH — I n)<*+ (1 j2 + /ab — I u)f3 

+ + I n)7 0-1) 
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On co vu lu ps for i und j lmve boon assigned, tin* 
expansion ovpr R involves assigning values to f u l 2 , 
ami Z 3 . Consider the Oth cell and let layers of 
cells be added to it layer by layer until it expands 
to fill all space. Let 77 designate a layer number 
starting at zero for the Oth cell only and taking 
values 1, 2, 3, and so forth as layers are added. 
The case for = 0 is a special case, as pointed out 
previously, and is computed separately. As 
values are assigned to 77, the following values are 
taken by l (there are three groups and all combina- 
tions in each group must be taken): 

(1 ) l x = v and —77 (2 values) 

1 2 — — r\ to 0 to 77, inclusive ( 2 rj-\-\ values) 

to 0 to 77, inclusive (277+ 1 values) 

(2) ^ = — (»| — 1 ) to 0 to (77— 1), inclusive 

(2?7 — 1 values) 

l 2 = 77 and — 77 (2 values) 

— 77 to 0 to 77, inclusive (277 4 1 values) 

(3) /, = — (77— 1) to 0 to (77— 1), inclusive 

(217—I values) 

/ 2 = — (tj — I ) to 0 to (77—1), inclusive 

(277— 1 values) 

/ 3 = 77 and —77 (2 values) 


Then 1 are 24 t / 2 + 2 possible combinations of /,, / 2 , 
and Z 3 . If there are r ions per unit cell, there will 
be t 2 — t terms in the total summation for the case 
of 77=0 and t 2 ( 24 t 7 2 + 2) terms for each case of 77^0. 

It is convenient to modify part of equation (lb) 
by changing from summation over R to summa- 
tion over 77. The following definition may be 
int induced: 

A=Z 2 1 (22) 

»j = 0 ij Clij 


In this equation the special provision for the 77 — 0 
case is tacitly included. A further definition may 
be given as 

a.=£ Z n- (2-i) 

7 j=o ij (l i j 

Then, as values are assigned to k, a value A« is 
computed. Tlie term A* will converge to A as 
k ^ 00. In practice, it is found that A* converges 
fairly rapidly and it is sufficient to let k go only 
as high as b or 8. A graphical extrapolation is 
then used to determine the limiting value of A. 
Once the limiting value of A is obtained, it is 


substituted into equation (lb) to give the coulom- 
bic energy as 

(24) 

The application of this method in computing 
the coulombic energy of the oxides of molybdenum 
is given in appendix B. 

POLARIZATION ENERGY 

The polarization energy is the energy resulting 
from the interaction of an electric field and a dipole 
induced in an ion or atom by that field. 

The polarizability a of an ion or atom is defined 
by the equation 

M=aF (2«) 

where 

fl electric- dipole moment 

- > 

F electric* field vector 

The energy is then given as (ref. 4 ) 

U p = — JaF-F (26) 

The electric field at a point in space due to a 
charged ion at another point in space is given by 
the following equation: 


F = 


Zel) 


(d-d) 


(27) 


where 


Z valence of the ion 

e electronic charge 

1) vector from the ion to the point of polariza- 
tion 


The energy for this euse is 


tt 1 72.2 d-d 

(2S) 

(d-d) 


or 


u 

j — >14 

(29) 

2 D 


Equation ( 29 ) shows the polarization energy to be 


proportional to the reciprocal of the fourth power 


G22837 — 62 


3 
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of tlu' distance from the point in question to a 
charge producing the electric field. Because of 
this, it is not considered necessary to account for 
elect rie-field-produeing ions at a great distance 
from t lii fc point of polarization. In computing 
the polarization energy of a “kind” of ion, a 
representative ion will he selected near the center 
of a cell. The electric-field vector is determined 
by taking a vector sum of the contributions of all 
ions in cases of 17 — 0 and r}— 1 in a manner some- 
what analogous to the coulombic-energy caleu- 
lat ion. 

Consider a polarizable ion designated by its 
position vector (q as follows: 


It — /,aa + / 2 b£ + /3C7 (d2) 

In the subsequent summation, values of / are 
assigned in the same manner as for the eoulombic 
energy and iti view of the short range of the 
polarize t ion effects only the cases for 17 — 0 and 
rj~ 1 are included. 

The vector P|j is found to be 

I) IJ =C1 J — T,— K (:w) 

This vector is associated with an electric-field 
cont rib 1 1 ion F u and 


( 1 j - gj 1 - a + g j 2 £+ g J3 7 0*0 ) 

This ion will be selected to represent one “kind” 

of ion designated by j. It is selected from all the 

- > 

ions 1, in a cell. The polarizing ions are the 

remaining ions Ij as shown in the following 
equat ion : 

Ii 2 jft+li 37 (d 1 ) 

The cell containing I, is positioned relative to the 
cell containing (q by the vector R as follows: 




0 * 4 ) 


1 vDjj-Djj/ 

If I)n is written in the form 


P] j — Piji«d"Pij20+Pij:i7 

(:«>) 

P iJl ~ Kn Pi 


Pij2— £12 1 12 ^>b r 

( 37 ) 

Pij3 — gj3 I J 



The expanded expression for Fj is 
Fj 


^ Z L (D| jiQ!-r^ij2^ ~l~I hjs7) 

° ^0 i (I)| J T + I), j r: + I)i^+2D 1 J 1 1) 1J2 cos 7 + 21)ij.d ) lj;} cos a+2D|j,Djj, cos 


(TS) 


There is one term in the above summation where 
— ♦ ► 

( f j = 1 1 and this must be discarded. If there are 
r ions per cell, then there are ( 27 t — 1) terms in 

the summation. The vector Fj comes out as a 
simple vector and from it and the polarizability of 
the polarizable ion, the energy of polarization is 
readily computed. The following modified expres- 
sion gives the polarization energy on a per mole 
basis for one “kind” of ion: 


where N Av is Avogadro’s number. There will be a 
polarization energy computed for each ion in a 
mole o* stoichiometric unit and these must be 
added to obtain the total molar polarization 
energy. 

In pnetiee, it is convenient to define a quant itv 
( A p )j a^ follows: 



( 40 ) 
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Tilt'll , 

i ? i=o(A p ) j (41) 

The. term (a;;)j muv ho obtained by modifying 1 > tj 
as follows: 

Du=- D,j (42) 

P 


Where P 7is a scaling parameter which may lx* 
taken as the shortest distance between two ions. 
Then, 


(a;:). 


(Dfj.Dr,)"* 


Th is treatment gives a more general character to 

the CO, expression and is somewhat analogous to 
the Madelung constant. For computing repulsion 
(MH*rgy suhse(juent ly it will he convenient to make 
the following definitions: 



The following energy expression is then written 
for one mole of one “kind” of ion: 


in appendix C. The inhomogeneity of the electric 
field was ignored because it is too complex a 
problem to bo included in this work. It is implied 
in reference 5 that the error introduced is not 
large. 

The application of this method in computing 
the polarization energy of the oxides of molyb- 
denum is given in appendix B. 

VAN l)ER WAALS ENERGY 

The Van der Wauls energy is a phenomenon 
of cohesion which has been explained by London 
(ref. (>) as the result of interaction of rapidly 
revolving dipoles in two atoms. The revolving 
dipoles are produced by the extranuclear electrons 
rapidly orbiting about the nucleus. When two 
atoms are close together, the revolving dipole in 
one synchronizes with the revolving dipole in the 
ot her in such a manner as to lower the total energy, 
and cohesion results. 

The mathematical development of the Van der 
Waals energy equation [‘('suits in a sum of de- 
creasing terms. Only two of these energy terms 
are considered herein: the dipole-dipole energy 

0<h (7)}, 



P 


and the dipole-quudrupole energy (eq. (S)), 


v„=2 (a,),— N (a p ), ■ (a p ) j 

j - j 

(45) 

or 

Up— (*«) 

This treatment of polarization involves some' 
approximations which should be pointed out. 
The energy of polarization calculated by the 
method just presented is based on the assumption 
that the electric field is homogeneous and that 
the polarizability of an ion is a constant and, in 
particular, not a function of the electric field. Tn 
the case of a real crystal, such [is those studied 
herein, the electric field is not homogeneous and 
tlu i polarizability depends, to some extent [it 
least , on the environment of the ion or atom. In 
this investigation, an effort was made to obtain 
polarizabilities of the ions for conditions similar 
to those existing in \Io() 2 and Mo() ;j , as discussed 



In these expressions, J and L an* constants for a 
particular crystal and p is the stalling parameter 
which is usually taken as the closest interatomic 
distance. In order to obtain the total molar Van 
der Waals energy, a summation must be taken 
of all nearest neighbor energy contributions for 
each ion in a stoichiometric unit. 

Consider an ion of ('lenient j and another ion of 
element k separated by a distance di. The Van 
der Waals energy expressions are 


(l dd)ijk — 

-Tjk 

<li" 

(47) 

(I d(j)l}k = 

w 

d, s 

(48) 


Then, summing over all neighbors k of ion j at 
distances <h and finally adding the contributions 
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of all “kinds” 
Waals energy 

of ions j gives t he molecular Van dor 
as 


TT V* ^ 

Udd— 2-1 A ti 

ijk Oi 

(49) 


T T v 1^* 

L dQ .1 8 

ijk <U 

(50) 

The scaling parameter p may be factored y it 
final expressions (eqs. (7) and (8)) 

hling 


u -=-? 


where 




li 

(r,i) 

and 

T . V 

h - 5 '<<>)* 

(52) 


<ir=-- <ii 

P 

(55) 


Th(» di polo-dipole coefficient, lias been derived 
by Kirkwood (ref. 7) as follows: 


may lx substituted for Kj and K k is determined 
in a similar manner with the following result: 


ii ( A) i \j 
(‘ 2 \ N .a j N k a k / J k 


The application of these equations in computing 
the Van dor Waals energy of the oxides of molyb- 
denum is given in appendix B (with the data, 
obtained being given in tables 1 and 2). 


REPULSION ENERGY 


The repulsion energy is the energy resulting 
from the displacement of two atoms or ions toward 
each other against an inherent force which tends 
to move them apart. This repelling force pre- 
sumably acts in accordance with the Pauli prin- 
ciple. There is no general fundamental expression 
that expresses the energy of repulsion between two 
atoms. There arc two empirical expressions in 
general use which are justified by the fact that 
they give good results. 

Consider t he case of a simple NaCI type crystal. 
The two repulsion energy expressions are, from 
the exponential law, 



a polarizability 

N number of polarizable electrons in the ion 

Equation (54) must be multiplied by e 2 a 0 1/2 to 
convert from atomic units to centimeter-gram- 
second units where e is the electronic charge and 
a 0 is the Bohr radius. 

The dipole-quadrupole coefficient has been ex- 
pressed in terms of t ho dipole-dipole coefficient by 
Huggins and Sakamoto (ref. 8, pp. 241-251) as 
follows: 



where 

X Av Avogadro’s number 
Kj ionic refraction 

fn order to express equation (55) in terms of 
polarizabilities, 


3R, 
4ttN Av 


(SO) 


and, from the inverse power law, 

( 59 ) 

In these equations B, g, and n are unknown 
parameter's. One of these, B, may be eliminated 
as follows: 

If the total energy is given by 

u=i T c +r P (fio) 

and it is assumed that, at d = R, 



where IL is the equilibrium interionic distance, then 
by substituting equations (58) and (59) into equa- 
tion (G( ) and applying the conditions expressed by 
equatio l (61), the energy expressions are 

u = u '( 1 -iTfi) m 

and 

1 r -(' •>■■■• 



THE CALCULATION OF LATTICE ENERGIES OF COMPLEX CRYSTALS 


19 


TABLE I.— VAN HER WAALS ENERGY DATA FOR Mo(L 


A 

Mo''* 


O, 


O n 


O,, 


Ion pair 


Bond 




Mam 

len ^th, 

O a H 

d ab 



Rah 



B 



1 


Or 

1 

1 . 04 



• 


Oi 

1 

1 . 08 



o, 

1 ,;, 

1 . 98 

1 . 23078 

0 . 1 1 382 

On 

!,; 

‘ 1 . 05 



o„ 

1 .. 

2 . 07 

■ 


On 

I 

l. 90 



Or 

1 

3 . 08 



Oi 

1 ■» 

3 . 04 



Oi 

1 

2 . 83 



Oi 

1 

2 . 74 

5 . 3207 1 

0 . 721935 

On 

E 

3 . 35 

i 


On 

1 

2 . 85 i 



On 

!•, 

2 . 77 


! 

Mo'' 

I ;> 

3 . 09 

\ 1 . 23078 

1 0 . 11382 

Mo iv 

I;, 

1 . 98 

] ] 

1 

<>i 

1 

1 . 49 



<>i 

1 <» 

2 . 91 



i Oi 

1/, 

2 . 45 



1 o, 

i. ( 

3 . 15 

► 5 . 32071 

► 0 . 721935 

On 

1 

2 . 85 



1 On 

1 

3 . 48 



<>n 


2 . i 2 



On 

l 

; 2 . 90 



Mo* 

H 

2 . 05 

l l. 23078 

1 0 . 1 1382 

Mo * 1 

l. ;> 

2 . 07 


1 


-S- >W'ab = 04:{9 „ ( y- MAjjdAB =() ,,,8859 

In Rab ab «»» 


These equations hold for the equilibrium distances 
only. The same approach that was used in the 
simple case is extended to the following complex 
case. The inverse power law is used because of 
its simplicity. 

Equations of the form of (59) and (63) are used 
as the basis of the computations of repulsion 
energy for the oxides of molybdenum because 
they are considered accurate at the equilibrium 
distance and because they are rather easily 
applied. The structure of these oxides is con- 
siderably more complex than that of XaCl and 
several modifications must be made. In the XaCl 
type crystal, only cation-anion contacts are con- 
sidered and the bond length is constant throughout 
the crystal. For the molybdenum oxides, both 


cation-anion contacts and anion-anion contacts 
are present. Also, the bond lengths range over 
several values for each type of contact. Pauling 
lias developed an equation to compute the repul- 
sion energy for the more complex crystals (ref. 2, 
p. 523). 

The potential energy between t he it h ion and the 
jth ion is given by Pauling as 


U 


u 


ZjZjc 2 

(In 


PftjB 0 e 2 



(64) 


The term on the extreme right expresses the repul- 
sion energy as 


F T r— 0ijB o e 3 


(r .+r,) D 

d n 


(65) 
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tabli: 2.— van m;u vvaals knkrcv data for. m 0 o. 


Ion 

pair 

M ab 

Bond 
length, 
B ah 

t 'AH 

( 1 a h 

A 

B 


Mo 

<h 

2 

i. 03 


' 



O, 

i 

Z 31 






On 

I 

2 . 08 

0 . 04050 


0 . 08451 1 


On 

1 

i. m 





O m 

1 

i. 82 

! 




<h 

Oi 

1 

2 . 58 

'' 





On 

2 

2. (>. r 






i »n 

2 

2. 8 ( 






On 

1 

2 . 81 


5 . 32071 


0 . 721035 


On 

1 

3 . 2 -t 






o,„ 

2 

2 . 40 







2 

3 . 40 





( h, 

o,„ 

1 1 

2 . 5 (‘ 






o„, 

1 

2 . 82 






* ( '"‘xr. 

! 

2 

3 . (H 


► 5.32071 


> 0 . 721035 



2 

2 . 77 






• - 


. 1 

... 



"Next layer. 

£ M ,V' Cab =0M()141 -*"-0.0061918 

iv ah \H 


where 

Bo a constant which is not dependent on the 
particular ion pair 
r i,r j ionic radii 

ii a constant which is dependent on the 
particular ion pair 

Also, 

/hj— l+pr+pr (06) 

where 

Z valence 

V number of electrons in the outermost shell 

for the particular ion 

In order to compute repulsion energy on a per 
molecule basis, a summation must be taken over 
all bonds associated with the molecule as 

W) 


where 

M u the number of bonds of type ij in a 

molecule 

d u = pi h 

The ccnstant B (> may bo eliminated in a manner 
similar to the NaCl case except that now con- 
sideration must be given to all kinds of binding 
energy present in the crystal. Thus (eq. (3)) f 


b t — l + T j, + l dd + l dq + l r 


where 


r A M e 2 Ca(‘ 2 J B , t> 2 V (c 

L t— ~ • , — — - S +B 0 e 2 2^ - u . (h; 

P P P P Ij P U,J 


G j j — Mu/ffjj 


(r i + r J ) n *j“ 


(09) 


Xow t ie conditions expressed by equation 
are applied as follows: 


dlh 

dp" 


\ M e 2 4Co , (‘ 2 6J sIj 


+ 


+ :r+5" B 0 e 2 2 . 

p P ij P 


n,,G 


m 1 ij 

»ii-n 


( 01 ) 

(70) 
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At equilibrium 



(71) 


TIip constant B„ may be factored from the re- 
sulting expression to give 


t> 1 / A M » 4Ca , 6 J , SL \ 

n , jG ij KZ t iW t («oV 


(72) 


?R„ n ‘i +1 


Then this equation may be substituted into tin* 
energy expression as follows: 




Gi 


TT — A M e 2 s y K 0 riiJ 

u * R~„ i ' , ‘ ^ ||„Q'„ 

ir Ro" i,+ ' 


a 


4Ca , 6.J 


<SL 


M | | "" I 

K „ 2 *' R,/’ c 2 RJ i ,2 R 9 




Cai ' 2 J L 

' r7 _ r7 - r ( ; 


In tliis rquat ion R„ has replaml p. 

By defining a repulsion coefficient p r ns 


(7:5) 


^ vr * (r.+r,) n «-‘ 
,)f - (r.+r,)^- 1 


(74) 


22 M u iii j /3ij 

ij 

equation (73) becomes 


r — A M e 2 Cac 2 J L , 

lt R 0 R 0 4 R»° R ,y [r 


/A M e 3 . 4Cae 2 , OJ .81 A , 7n 

(■«„■+ K/ + ic + ie) (,5) 


Thus, equation (4) is obtained as 

T T P = — p r t U c + 4 Up + 6 1 J<w + 8 1 r { i fl ) 

In computing the repulsion coefficient, it is 
necessary to assign a value to iii j. It has been 
found that n may be assigned a value of 9 and 
for most crystals any error introduced will be 
small. Pauling has determined a set of values of 
ii for ions that are isoeleet ionic with the inert 
gases (ref. 2, p. 509). For a beteronuelear bond, 
an average is taken to arrive at a value for lUj. 

The repulsion coefficient is computed lor the 
two molybdenum oxides in appendix B and this 
in turn is used to compute the repulsion energy 
(tables 3 and 4). 


SUMMARY OF ENERGIES 

The contributions to the lattice energy which 
have been calculated in appendix B are presented 
as follows and are compared with the experi- 
mental lattice energy determined by means of the 


Born- Haber cycle. All 

energy units 

; are kilo- 

calories per mole. 

MoOj 

MoG» 

Coulombic, i c 

—2, no. oo 

-5, 850. 00 

Polarization, U p _ 

-53. 12 

— 868. 46 

Dipole-dipole, U d(1 

D i pole-qu a< l rupole , 

— 105. 99 

— 72. 70 

l ( |q 

— 7b. 57 

— 25. 13 

Repulsion, U r -_ 

( calculated lattice 

402. 44 

1, 201. 12 

energy 

— J , 883 

-5, 555 

Experimental lattice 



energv 

Difference in calcu- 

— 2, 805 

— 0, 009 

lated and experi- 
mental lattice 



energv 

Ratio of calculated 

— 922 

-454 

and experimental 
lattice energy 

0. 071 

0. 924 


INTERLAYER ENERGY 

The bonding energy that exists between layers 
of a crystal is important as a factor affecting the 
strength and stability of the crystal. A good 
deal of insight as to the nature ol the interlayer 
bond may be obtained by computing the coulombic 
energy of an isolated layer and comparing it 
with the coulombic energy of the crystal. 

The layer coulombic energy may be computed 
in the same way that the crystal coulombic 
energy was computed except that the system 
for assigning values to l u / 2 , and / 3 must be 
changed and position vectors of atoms within a 
single layer only must be selected. In the case of 
molybdenum trioxide the layers are oriented 
parallel to the b-e plane. If position vectors are 
selected to represent atoms within a single layer 
in the crystal cell and the crystal is expanded by 
adding cells in the b-c plane only, then a coulombic 
energy may be computed by using the method 
described previously and this will be the coulombic 
energy of an infinite single layer. 

The following system of assigning values to 
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TABLE 3.— REPULSION ENERGY DATA FOR MoO* 


[For Mo**-: r--0.<>8, Z-+4, P=IO; for O 2 *: r= 1.40, Z -----2, P-8. Galculatod ]>,=(). 1295J 


Ion pair 

n i j 

Mii 

du 

i’i ~h r j 


i 

j 



Mo"' 

<>1 

8. 5 

1 

. 

I. 939 

i 

2. 08 

1 15 


<>1 

8. 5 

1 

1. 984 

2. 08 

1. 15 


<>1 

8. 5 

l / 

1. 982 

2. 08 

1. 15 


<>ir 

8. 5 

1/ 

■>. 054 

2. 08 

1. 15 


<>n 

8. 5 


2 . 073 

2. 08 

1. 15 


0,1 

8. 5 

1 

!. 992 

2. 08 

1. 15 

Oi 

Oi 

t 

1 ;, 

3. 081 

■2. 80 

0. 50 


o, 


l „ 

3. 939 

2. 80 

. 50 


O, 

i 

i 

2. 830 

2. 80 

. 50 


O r 

7 

1 

2. 742 

2. 80 

. 50 


On 

7 

H 

3. 347 

2. 80 

. 50 


Ou 

7 

1 

2. 849 

2. 80 

. 50 


On 

7 

D 

2. 705 

2. 80 

. 50 ! 


Mo" 

8. 5 

1 4 

3 089 

2. 08 

1. 15 


Mo iv 

8. 5 

h 

‘ l. 982 

2. 08 

1. 15 j 

! On 

Oi 

7 

1 

1 1. 480 

2. 80 

0. 50 


Oi 

i 

1 t 

2.910 

2. 80 

. 50 


Oi 

7 

i 

2. 440 

2. 80 

. 50 


<)[ 

7 


3. 140 

2. 80 

. 50 


On 

7 

1 

2. 845 

2. 80 

. 50 


On 

7 

*4 : 

3. 475 i 

2. 80 

. 50 


On 

i 

1 / ! 

2 . 7 1 7 

2. 80 

. 50 


On 

i 

1 

2. 890 

2 80 

. 50 


Mu' 

8. 5 

Lj 1 

2. 054 

2. 08 

1. 15 


Mo' 1 ' 

8. 5 


2. 074 

2. 08 

"■ j 


/i, I 2 , and A { is applicable to layers lying parallel 
to the b-c plane. The value of zero is assigned to 
(\ throughout the computation. The 77— 0 case 
is computed separately. For all other values 
assigned to 77, l must take the values given as 
follows (all combinations of each group must be 
taken) : 

(1) Ji — 0 

U~r) and — 7/ 

4= — 77 to 0 to 77 , inclusive 

( 2 ) /, = () 

1-2— ~(y~ 1) to 0 to (77— 1 ), inclusive 

/ 3= — 77 to 0 to 77 

In this computation, l lie following position vectors 
were taken: Mo H1 , Mo IV , Oj, 0\\ Oft, Oh, O/ft, and 
Ojh. Examination of figures 8 and 9 shows that 


these atoms lie within a single layer and single cell. 

The computation just presented has been 
carried out and the layer coulombic energy is 
given a ong with the crystal coulombic energy for 
compar son as follows: 

Laver U t .= — 5,741 kcal/mole 
Crystal U c = — 5,850 kcal/mole 

One ,vay that tlu* layer-structure energy rela- 
tionships may be interpreted is to compare the 
stabilization effect of neighboring atoms on a 
molecular unit. For t his purpose it would 
perlmp? be better to use a molecular unit in the 
eonfigui ation tliat it would assume isolated in 
space. This type of information is frequently not 
available. Therefore, a molecular unit will be 
used in the configuration that the unit assumes in 
a crystal. For these computations the atoms 
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TABLE 4.— REPULSION ENERGY DATA FOR MoO, 


l For Mo' 11 : r -04)2, P = 8; for O 2 -: r-1.40, / m 2, P = 8. Calculated p r - 0. 1 2(i(>] 


Ion pair 

i i 

n.j 



Ti + Tj 

flu 

Mo 

Oi 

8. 5 

2 

t. 94 

2. 02 

1. 50 


Oi 

8. 5 

1 

2. 4 1 

2. 02 

1. 50 


On 

8. 5 

1 

2. 08 

2. 02 

1. 50 


On 

8. 5 

1 

1. 89 

2. 02 

1. 50 


Oin 

8. 5 

l 

l. 84 

2. 02 

1. 50 

<>i 

o, 

/ 

l 

2. 58 

2. 80 

0. 50 


On 

7 

2 

2. 05 

2. 80 

. 50 


On 

- 

2 

2. 80 

2. 80 

. 50 


Or, 

7 

l 

2. 81 

2. 80 

. 50 


< >ii 

7 1 

1 

4. 24 

2. 80 

. 50 


Oiri 

t 

2 

2. 49 

2. 80 

. 50 


On, 

i 

2 

4. 49 

2. 80 ; 

. 50 

On 

0 1 1 1 

t 

I 

2. 50 

2. 80 

0. 50 


On, 

i 

1 

2. 82 

2. 80 

. 50 


On, 

7 

2 

4. 00 

2. 80 

. 50 

Dm 

Onr 

7 

■ 

2 

2. 77 

2. 80 

0. 50 

L 



L 

1 

... _ ... 

. ... .. j 


Mo 111 , Oi, On, one I Oi'ii wore selected and the 
electrostatic potential energy of this group was 
computed by using the position vectors to estab- 
lish t he configuration. The electrostatic potential 
energy or coulombie energy of this group was 
found to he —4,984 kcal/mole. 

A stabilization factor may be defined as 


layer stabilization factor is 98.1 percent of t lie 
crystal stabilization factor. 

A different way of interpreting the layer 
coulombie energy is to compare differences in 
energy rather than ratios of energy. Stabilization 
energy may be defined as 

Crystal U„=Crystal r c — Molecular V c (77) 


Stabilization factor = 


Crystal U c 
Molecular V v 


(7b) 


and similarly defined for the layer coulombie 
energy. This stabilization factor represents the 
increase in the magnitude of coulombie energy 
gained by a molecular unit when it is located in an 
infinite crystal rather than when it is isolated in 
space. The stabilization factors for the crystal 
and layer are given as follows: 

(Tvstal stabilization faetor= ^^^ — 1.174 

— 4,984 


Layer stabilization factor— 


-5,741 

-4,984 


.152 


On comparing these values, it can be seen that the 


and similarly defined for the layer stabilization 
energy. Values obtained with this method are 
Layer l s — — 757 kcal/mole 
Crystal lJ a =— 806 kcal/mole 

From these values it can be observed that the 
layer contains 87.4 percent of the crystal stabili- 
zation energy. 

The difference between the crystal coulombie 
energy and the layer coulombie energy is —109 
kcal/mole. This value represents the coulombie 
energy that a molecular unit gains when an infinite 
crystal is placed (with the right orientation) 
against each surface of the layer. Since two 
surfaces are involved, half the difference between 
the energies may be considered to represent an 
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average interlayer coulomhic energy. 'Finis, the 
average interlayer eoulombir energy is —54.5 
keal/mole. 

From table 2, it may be seen that there are two 
Oni — Out and two On — Om next-neighbor type 
interlayer bonds per molecular unit. From these 
data, the interlayer Van der Waals energy may 
be computed and the values obtained are given 
with the average interlayer coulomhic energy for 
comparison as follows: 

Average interlayer U c =— 54.5 keal/mole 

Average interlayer U dd =— 9.2 keal/mole 

Average interlayer l dq = — 2.0 keal/mole 

These values indicate that t li o interlayer bonding 
mechanism is not predominantly Van der Waals 
type but that the bonding energy is small never- 
theless. 

COVALENT EFFECTS 

In reviewing the summary of energies, it is 
readily apparent that the electrostatic model of a 
crystal is capable of accounting for the total lattice 
energy of molybdenum trioxide although it is 
entirely possible that there is cancellation of errors 
in summing the different- energy contributions. 
It is also apparent tlmt this model does not ac- 
count for t lie total lattice energy of molybdenum 
dioxide. A quantitative estimate of covalent 
energy contributions will not be undertaken here- 
in but it is worthwhile to examine possible evi- 
dence of covalent effects. Consider first the 
dioxide* of molybdenum. According to Pauling 
(ref. 2. p. 4^7) the two d-electrons on the Mo*" 
ion interact between pairs of molybdenum ions 
and form a bond of bond number 1.47. 'Flu* alter- 
nate close and far spacing of the molybdenum ions 
in Mo() 2 plus the (dose bond distance of 2.49 A 
as compared with the metallic bond distance of 
2.72 A tends to support this hypothesis. It is 
difficult to estimate accurately an energy contri- 
bution from this source. The energy of atomiza- 
tion of molybdenum is 155.5 keal/mole. If 
molybdenum is considered to have a metallic 
valence of 0, then then 1 are three single bonds per 
atom and the single bond energy would be 51.8 
kcal/bond mole. In the \lo() 3 crystal there is 
one half-bond of bond number 1.47 per mole. 
A reasonable estimate of an upper limit to this 
energy contribution would be 100 keal/mole and 
this value is far short of the 922 keal/mole of 


energy which is as yet unaccounted for. Other 
bond distances in Mo() 2 may be examined and 
compared with the following covalent bond 
distances: 

From the Shomaker-Stevenson relationship dis- 
cussed in reference 9, 

Bond Mo — (), 188 A 


and from the 0 2 molecule, reference 10, 
Bond ()-(), 1.207 A 


The ionic bond distances (based on 0.08 A as 
Mo H radius) are 2.08 A and 2.80 A for Mo~() and 
O—O, respectively. In the octahedron, the Mo— O 
bond distances are as follows: 1.94, 1.98, 1.98, 
1.99, 2 00, and 2.08 A. The 0 — 0 distances vary 
from 2.08 to 4.00 A. These distances do not give 
a conclusive indication of the covalent or ionic 
nature of the bonds, especially when allowance is 
made for distortion of the lattice and polarization. 

A possible mechanism which could account for 
the energy in Mo() 2 would he an electrostatic 
energy contribution induced by a covalent bond 
bet wee i the two close molybdenum ions. If a 
covalet t bond is present between these two ions as 
Pauling suggests, then the electron density will he 
increased bet. ween the two ions and decreased at 
the encs of this ion pair as compared with the case 
of two individual ions. A crude model was set 
up by placing a charge — q midway between every 
close pair of molybdenum ions in the Mo() 2 lattice 
and a c large was added at the position of every 
molybe enum ion. The quantity A from equation 
( 22 ) 


A= 


=£ s v f 

ij 0(j 


and fix in 


Ur 


195.9 A 

£ 


(78) 


was computed for the cases of <[ = 1 , 2, and 4. 
The results of this work are shown in figure 10. 

On tlie basis of coulomhic energy alone, with 
1 00 kcr 1/mole being allowed for the Mo — Mo bond, 
a value of A — 70.7 is needed to account for all 
the energy in the Mo0 2 lattice. This value of A 
is equiv alent to q — 0.4. From a qualitative con- 
cept of a covalent bond of bond number 1.47 
and in consideration of the crude nature of this 
model, this seems to be a reasonable value for q. 
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220 r DISCUSSION 



Summarizing the consideration of the Moth lattice 
shows that the evidence for a covalent bond 
between the close pairs of molybdenum ions is 
good. There is probably some covalent contri- 
bution to the energies of the Mo — 0 bonds but 
there is no conclusive evidence that this contri- 
bution is large. 

In the case of the trioxide of molybdenum, 
there is no evidence, other than possibly a bond 
distance, of an appreciable covalent contribution 
to the bond energies. A comparison of bond 
distances shows some Mo — () bonds to be at the 
right distance to be covalent. However, the 
large amount of polarization present in tins 
lattice makes it possible 1 for some ionic bonds to 
exist at distances other than the sum of ionic 
radii. This possibility tends to reduce the impli- 
cation that these bonds must be considered 
covalent on the basis of distance alone. 


MODEL 

The model used in these calculations has 
already been described; however, there are 
several comments concerning the use of this 
model as applied to the more complex crystals 
studied hero in contrast to the simple crystals 
studied prior to this investigation. 

The idea that the net charge of an ion can be 
considered to be located at a point in the center 
of tire ion comes from the approximation that 
the electronic charge in an ion has a spherically 
symmetrical distribution about the nucleus either 
in a series of thin shells or a solid distribution 
with perhaps some radial variation in density. 
It can be shown in classical electrostatic compu- 
tations that the electric field generated by this 
model of an ion will be the same as that generated 
by a point charge located at the center of the ion. 
The requisite key to this concept lies in the 
spherical symmetry of the charge distribution, 
and the electric field resulting from the net ionic 
charge is considered only at points outside the 
ion. 

The electrostatic model has been applied to 
many crystals with great accuracy. The crystals 
that have been amenable to methods of calculation 
based on this model have a number of unique 
features in common. The geometry of the 
crystals was simple and highly symmetric. The 
ions wen 1 isoelect ronic with the inert gases. 
The positive ions were highly electropositive and 
the negative ions, highly electronegative. The 
polarizabilities of the ions wore relatively low. 
Although these features appear to be desirable 
for accurate use of the electrostatic model, it 
is apparent that there can be appreciable deviation 
from these features and still obtain calculated 
energies that are within 90 percent of the ex- 
perimental energies. This value may be com- 
pared with about 98 percent for the more ideal 
crvst a Is. 

If the concepts of electron charge distribution 
that- are derived from quantum chemistry are 
representative- of a true physical situation, then 
the real atom does not possess the kind of charge 
distribution needed in the electrostatic model in 
an exact sense, although the approximation must 
be very good. The features of the ideal ionic 
crystal appear to be such that they accomplish 
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one basic act. They tend to preserve the sym- 
metry of charge distribution of an ion about the 
nucleus. The large difference in electronegativity 
of the cations and anions assures complete charge 
transfer leaving the ions isoeleetronie with inert 
gases. The symmetry of the crystal provides 
that the electric field at the position of each ion 
is zero and the distribution of charge in an ion 
is then not altered by polarization. 

The application of this model to a distorted 
crystal like molybdenum trioxide incurs several 
effects worth observation. In this work it was 
assumed that the position vector of an ion es- 
tablished the position of the nucleus of the ion 
and that this was the center of net charge. If 
this assumption is true the electric field ■which 
has been shown to exist at the nuclei would cause 
a polarization that would destroy the symmetry 
of the ion charge distribution. This polarization 
would induce two effects: (1) there would be 
established a polarization energy such as that 
computed in this work and (2) the symmetry of 
charge distribution would be changed and the 
center of positive charge would no longer coincide 
with the center of negative charge. In order 
for the model used in this investigation to ap- 
proach the real physical situation more exactly, 
each ion should he considered not as a point 
charge hut rather as two point charges repre- 
senting the centers of positive and negative 
charge. This representation is still an approxi- 
mation and is perhaps a less accurate approxi- 
mation to a real polarized ion than the point 
(diarge approximation to a nonpolarized ion. It 
Is probably true that the separation of centers of 
charge is very small because the inner electrons 
of an ion are so tightly bound although the orbitals 
representing the outer electrons may be quite 
distorted. The position vectors probably do 
locate the positions of the nuclei with good 
accuracy, as far as the model is concerned, since 
the position vectors are established by X-ray 
diffraction maxim urns which are a function of 
electron density. Even for polarized ions, the 

electron density should be greatest near the 
nucleus. The distortion of the outer orbitals 
should have some influence in the measurement 
of atom positions in a crystal and this is possibly 
an important factor in barring great accuracy 
in this type of work. 

A distortion of t ho outer orbitals of an ion w ould 


be expected to have an effect on the computation 
of Van der Waals energy and repulsion energy. 
The polarizability of a distorted ion would be 
expected to be nonisotropic. Also, the compressi- 
bility of a distorted ion would be expected to be 
nonisotropic. The Van dor Waals energy and 
the repulsion energy of individual bonds are 
probably not us accurately computed in these 
crystal:, as they would have been had the crystals 
been symmetric and without- distortion of the 
ions. Since there are bonds in all directions about 
a distorted ion, it is perhaps reasonable to assume 
t fiat the effects of this nonisotropy on calculations 
involving the whole crystal will cancel in an 
averaging process. 

In many studies of bond types, much emphasis 
is placed on the bond length or distance. Atoms 
and iors have been assigned various radii, that is, 
eovalei t, ionic, metallic, Van der Waals, and so 
forth, md these are used with the assumption 
that a bond distance should be equal to the sum 
of the appropriate type of radii. In this sense 
distance establishes bond type. These radii are 
determined from measurements made on crystals 
or moli rules known to possess certain types of 
bonds. The distance relationships appear to 
work o it very well on many substances and the 
bond t} pes determined by this method correspond 
qualitatively with other features of the substance. 
It is tc be expected that these relationships will 
not be valid for substances containing distorted 
ions all hough, with some allowance for the dis- 
tortion, bond distances may be compared with 
sums of radii to provide a rough indication of bond 
type. 

When the types of bonds in a crystal are being 
considered, comparisons with radii sums should 
be mat e with some allowance for the direction 
and nn gnitude of the electric field at each ion, 
the polarizability of the ion, and the expected 
shape < f the ion. An analysis of bond distances 
is made in t he section entitled “Covalent Effects M 
where he bond distances in Mo() 2 are used to 
obtain some insight as to bond type. On the basis 
of polarization energy, \lo() 2 is considerably less 
distorted than Mo() 3 and the bond-distance 
criteria is useful to some extent. The bond 
distances in MoO s do not appear to be very 
meaningful as to bond type. Of course, in con- 
sidering bond distances in crystals, it is necessary 
to be certain that the geometry of the crystal 
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permits ion contact at the point that the bond 
is under consideration. 

EXPERIMENTAL WORK 

Since the energy computations are- based on 
the relative positions of the atoms or ions and 
these positions have been experimentally estab- 
lished, the importance of the accuracy of these 
measurements is worth consideration. 

In X-ray diffraction work, the size and shape 
of the crystal cell may be determined with great 
accuracy, that is, five or six significant figures. 
The positions of the atoms in lire cell are deter- 
mined with less accuracy, perhaps three to five 
significant figures. The data used in t his investi- 
gation wore obtained by a highly experienced 
group of specialists using single-crystal methods. 
In general, the positions of the molybdenum 
atoms are believed to be established with accuracy 
at least- as good as three significant figures. Tire 
positions of the oxygen atoms arc* considered to be* 
established with only fair accuracy. Although 
three significant figures wore assumed in estab- 
lishing the positions of the oxygen atoms, the data 
in the original reports (refs, tl and 12) imply that 
this is perhaps a bit optimistic for some of the 
oxygen positions. 

\o effort was made to determine the effect of 
atom position accuracy on the energy computa- 
tions. The displacement of an ion from its true 
equilibrium position in a crystal would have the 
effect of increasing the lattice energy. Tin* extent 
of the increase would depend upon the displace- 
ment-, or error in position, and the shape of the 
potential energy “well.” It is not necessarily 
true, however, that a real crystal exists in an 
idealized state of lowest potential energy espe- 
cially at moderately high temperatures and it is 
possible that the displacement of an ion could 
lower tin* lattice energy. It is not known what 
effect an error in an atom position vector would 
have on the lattice energy. This effect could lx* 
determined but it would involve much computa- 
tional work. The experimental work on which 
these computations wen* based has been accepted 
in this investigation as being accurate. 

TEMPERATURE 

There are temperature effects present in the 
energies discussed herein which have been ignored 
in this work hut are sometimes considered in 


other studies of this kind. The lattice energy of a 
crystal should properly he the electrostatic poten- 
tial energy of the ions in a crystal at 0° K. In 
some work, even the zero-point vibrational energy 
is removed. The error involved in ignoring room- 
temperature effects is of the order of magnitude of 
5 keal/mole. 

The Born- Haber cycle calculations an* effec- 
tively isothermal at 298.15° K. It would be 
possible* to determine the measured lattice energies 
of 0° K but there are other considerations which 
make this hardly worthwhile. 

The calculated lattice energy determined in 
this work is tlu* electrostatic potential energy of 
the ions in a crystal without any temperature- 
induced kinetic energy effects. Even tin* zero- 
point vibrational energy has been omitted. The 
crystal lattice* is not that winch would naturally 
occur at 0° K but it is a lattice whose configuration 
is established at, presumably, room temperature' 
and the atom positions are mean positions. 

Those temperature* corrections would appear in 
this work in comparing the calculated lattice 
energy with the* measured lattice energy. It was 
felt that the magnitude of tl ic correction was 
trivial; the effect is compensated somewhat by 
the omission of the* AnRT term which was also 
ignored in the Born-1 laber cycle calculations. 

CONCLUDING REMARKS 

A method has been developed which makes pos- 
sible the calculation of the components of the 
lattice energy of an ionic crystal; this method has 
been applied to molybdenum dioxide and molyb- 
denum trioxide. The method is unique in that 
it may be applied to complex and distorted 
crystals. It may also be used to study perturba- 
tions in a crystal lattice. Essentially all that is 
necessary for this method to bo applied is that 
the crystal lx* ionic and possess periodicity from 
cell to cell. The main contribution from this 
research is the development of means of comput- 
ing coulombic energy and polarization energy. In 
addition, established techniques for computing 
Van der Waals energy and repulsion energy have 
been modified and extended to provide a complete 
method for studying ionic crystals. 

This method was applied to two oxides of 
molybdenum which are highly distorted struc- 
tures. The trioxide of molybdenum was found 
to he amenable to this method of study. The 
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calculated lattice energy was 92.4 percent of the 
measured lattice energy. Molybdenum trioxide 
is known to have a layer structure and, as would 
be expected for this type of structure, it was found 
to have a relatively high polarization energy. The 
interlayer bonding energy was studied and the 
interlayer bonding was found to be relatively 
weak. This weakness would account for the high 
vapor pressure of molybdenum trioxide. Although 
the interlayer bonding is weak, Van der Waals 
bonding is not predominant as might be expected; 
rather there is a eoulombie contribution which 
predominates. 

The dioxide of .molybdenum was found not to 
be amenable to this method of study. Only 07.1 
percent of the measured lattice energy was given 
by the energy calculations. A metal-to-metal 
bond involving d-electrons in pairs of molybdenum 
atoms was investigated by a simple model which 
contains a hypothesized induced electrostatic 
effect. This model is very crude but the calcula- 
tions imply that the metal-to-metal bond exists, 


and the induced electrostatic effect- will account 
for the energy discrepancy. 

The relationships between energy-producing 
mechanisms and the crystal structure of these 
oxides appear to be feasible according to general 
concepts of structural chemistry. Much doubt 
has existed as to the nature of the bonds in oxides 
like molybdenum trioxide. Some 1 investigators 
have expressed the opinion that these compounds 
are covalent. This research indicates molyb- 
denum trioxide to be an ionic compound. The 
nature 1 of the bonds in the dioxide of molybdenum 
is not clearly established in this investigation. 
The evidence is moderately strong for ionic bonds 
plus the special atom-pair metal bond. It may 
be speculated that many of the* distorted metal 
oxides will be found to be of an ionic nature with 
perhaps additional anomalous mechanisms respon- 
sible for their unique character. 

Lancle v Research Center, 

National Aeronautics and Space Administration. 

La> iu.ky Air Force Rase, Va., March i\ f ( Jti / . 



APPENDIX A 


CRYSTAL-STRUCTURE CALCULATIONS 


For a number of years X-ray diffraction tech- 
niques have been employed by specialists to deter- 
mine the structure of crystals. By means of these 
measurements the type of crystal, the space group, 
cell dimensions, atom position parameters, and 
other feature's may be determined. With modern 
equipment and techniques these measurements 
may be made' with great accuracy especially where 
single-crystal measurements are made. In the 
following section on structural data the experi- 
mental crystallographic elata on molybdenum di- 
oxide' and molybdenum trioxide are given. These 
data are 1 taken from the original papers and coin- 
pendiums on the 1 structure of crystals. The sec- 
tion entitled “Atom Positions and Bonds” gives 
computed data in the form of position vectors, 
bond vectors, bone! length, and bond angle's. This 
is the detailed numerical information derived from 
the X-ray diffraction measurements which estab- 
lishes the structure of these crystals, implies the 
nature of the chemical bonds, and forms the basis 
of the energy computations. 


CRYSTAL-STRUCTURE DATA 


Molybdenum dioxide: 

The X-ray diffraction data for molybdenum 
dioxide (ref. 11) are presented as follows. The 


notation is that used in 
( Yvstal type 
Space group _ 

Coll content 
Cell dimensions: 

o 

H. A 

o 

I), A 

o 

e, A _ 

0, deg o 

Cell volume, A : f 

Density, g/cnP 

Coordinates of equivalent 

4(e) positions. 


'forence 4. 

Monoclinic 

n-P2, 

4 units MoOo 

0.584 

4.842 

5.fH)8 

120.04 

120.4 

7.02 

posit ions: 

x, y, z 
x, y, z 
x, fs + v, U — z 
r \, ! 2 — v , U+z 


The atom position parameters are: 

For four Mo in 4(e), 

x — 0.242, y= 0.000, z- 0.017 

for four ()i in 4(e), 

x — 0 . 1 I , y = 0.2 1 , z — 0.24 

and for four On in 4(e), 

x — 0.49, y— 0.70, z = 0.40 

Molybdenum t rioxide: 

The X-ray diflVaetion data for molybdenum 
trioxide (ref. 12) are presented as follows: 

Crystal type .Orthorhombic 

Space group DJJ— l’nmu 

( Yll content _ _ . 4 units MoO s 

( Yll dimensions: 

a, A . 

b, A - 

c, A . 

( .''ell volume, A : V 

Density, g/cnP 

Coordinates of equivalent positions: 


14. 85 
4. 090 
4. 900 
204. 0 
4. 50 


4(c) positions. 


X Y z 
C-x, )H'Z 
'■i + x, Ya> C-z 


The atom position parameters are: 
For four Mo in 4(e), 

x = 0.0998, z = 0.084 

for four Oj in 4(e), 

x = 0.445, z = 0.525 

for four On in 4(e), 

x=0. 100, z = 0.50 

and for four () m in 4(e), 

x — 0.240, z-0.015 
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ATOM POSITIONS ANI) BONDS 

The X-ray crystallographer presents structural 
data in the concise form given in the previous 
section. These data have been converted to 
a more usable form by the method described in 
the section entitled “Crystal-Structure Calcula- 
tion.” By tliis method the position of each 
atom in a crystal can be expressed by a position 
> * 
vector 1. The bond between two atoms and 

lj may be expressed by a vector as follows: 

idVii-c 

The bond distances L are determined by the dot- 
product relationship 

The interbond angles are also determined by the 
dot-product relationship; for example, in the 

case of ldj and 1 d+ 

Zl jl |Ij=0jU 

..(iXui) 

= <•‘08 — T 

'Hie notation used is keyed to the figures and 
the particular atoms may be located by referring 
to the vectors identified by element symbol, 
subscripts, and superscripts. 

Molybdenum dioxide: 

The atom position vectors for molybdenum 
dioxide have been computed as an example to 
illustrate the method in the section entitled 
“Crystal-Structure Calculations.” There are four 
equivalent positions, one “kind” of molybdenum 
atom, and two “kinds” of oxygen atom. The 
reference system is nonorthogonal. The atom 
position vectors are as follows: 

Mo»~* = 1 ,295a + 0.000/3+ 0.09537 
Mo^ = - 1 .295a — 0.000^3 — 0.095:17 
Mo^= - ] .295a + 2.42 10+2.7O9 t 
M o 1 ^— 1 .295a + 2.42 10+ 2.8997 


01 = 0.51 4a + 1 .0170+ 1 .8407 
0!'= —0.614a— 1 .0170- 1 .840 y 
(); ii =-o.Gi4a + :k4:kS0~f 1.4587 
(If = 0.0l4a+l. 4040+4 I0O7 

Oh = 2. 178a + 8.8890+ 1 .0827 

Oh- -2.1 78a - 8.8890 - 1 . 0827 

()fi - - 2 . 1 7 8a + 5 . 8 1 ( )0 + 1 . 122 7 

Oh = 2. 178a -0.9080 + 4.4807 

These position vectors have been used to sketch 
projections of the lattice. The projections are 
shown in figures 0 and 7, It may be observed 
that if i lie c-axis in figure 0 is drawn perpendicular 
to the a -axis, there will be only slight deviation 
from true periodicity due to the distortion of the 
lattice and from this viewpoint the lattice is very 
similar to the rutile lattice. 

The bond vectors have been computed and are 
listed with the bond distances. These vectors 
have been computed in an orthogonal reference 
system for convenience. The transformation 
used is 

* 

a — l 

A. 

0 = ] 

7= -(cos 59.06°)!+ (sin 59.()G°)k 

The transformed vectors may be used as long as 
the bounds of a single crystal (‘ell are not exceeded. 
Periodicity cannot be carried over in the trans- 
format on. Along the a-axis the alternating 
Mo-Mo distances are 2.497 A and 8.090 A. 
These distances are important in postulating a 
unique Mo -Mo bond. With appropriate vector 
subtraition, vectors to represent bonds are ob- 
tained as follows: 

Mo“iJi v =1.17i-1.02J + 1.24k (L=1.9S A) 

Mo^of = 1 .82! + 1 .02] - 1 .07 k (L=1.9S A) 

M+ 1* = — 0. 1 1 St + 1 .40] + 1 .88k (L - 1 .94 A) 
Mo“*3h~ — 1. 58! + 0.96Sj— 0.881k (L=2.06 A) 

]VXo“i3fc= — 1.511 — 0.96SJ+ 1.045k (L=2.08 A) 
-O.OOSi- 1 .45] - 1 .80k (L= 1 .99 A) 
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( )y*()f= -o. i 5t>i-2.<>34) +2.3 1 k 
( )"< )f = 1 .291 - 2.42 j — ().()9<>k 
Oi.of = 2.751- l .98J + 2. 12k 
( )|,( V = 2.<>8i — 0.049] - 0.1 92 k 
()ri()? = 1 .241 + 0.43(i) + 2.(>0k 
( )('()+- 1 .441— 0.387} — 2.4 1 k 
( )i,()f = 2.9 11 + 0.049) — 0.1 92k 
( )i ; ,( )? = 2.831 -)• 1.99J— 2.12k 
( );■;( )+= 1.391 +2.47) + 2.8sk 
<)i,()y= 1.471 + 0.430) + 2.2 lk 
)y 1 .391 4-2.37) + 0.288k 
OyjOy = — 0.0501 +2.8(1) + 2.69k 
()y,0^= -0.0751+1 .94 j — 1 .93 k 
011(^1= - 1 .521 + 2.42) + 0.480k 
O'/jOy, = - 1 .441 + 0.4S5] 4- 2.4 1 k 


(L = 3.082 A) 
(L— 2.74 A) 
(L=4.00 A) 
(Ij— 2. OS A) 
(L=2.91 A) 
(L = 2.82 A) 
(I..=2.92 A) 
(L=4.0« A) 
(Ij=2.S5 A) 
(Ij=2.09 A) 
CL =2. 77 A) 
(L=3.93 A) 
(L=2.73 A) 
(L=2.90 A) 
(L=2.S5 A) 


The bond length in each case is shown on (he 
right. The bond vectors may point toward either 
end of the bond blit the direction must, be con- 
sidered in computing the interbond angles. The 
bond vectors given here point from the atom on 
the left toward the atom on the right. 

The interbond angles are computed using the 
dot-product relationship. The angles with apexes 
at the AIo lu atom have been computed and arc 
tabulated as follows: 


ZOyMt+Oy =88.70° 
ZOiMo iii Oi,=N4.<>5° 
Z()i l Mo i “()“i=91.3S° 
ZOyj I Mo"'Oi' r = 94 . 1 7 ° 
ZOi+Yh+Oy =86.9(1° 
Z()yMo i “()y i =92.25° 
ZGyMi+Oyi =91.(55° 
ZOy; Mo"‘Oyi = ss.s2° 
ZOy I Mo i, ‘Oy v =S2.81 ° 

zoy v Mo iU oy i = 92 . 0 (> o 

Z()y i Mo iii Oi I =92.19° 
ZOi I Mo iii Oy i = S2.87° 


Molybdenum trioxide: 

The atom position vectors for molybdenum 
trioxide are determined by using the method 
described for molybdenum dioxide. There ate 
four equivalent positions, one “kind” of molyb- 
denum atom, and three “kinds” of oxygen atoms. 
The reference system is orthogonal. The atom 
position vectors are as follows: 

Mo* = 1 .381 ‘ 0.924) + 0.333k 
Mo“== - 1 .381 + 2.77) — 0.333k 
Mo a, = 5. 541 + 2. 77) + 2.32 k 
M, + = 8.301 + 0.924)+ I ,05k 

Ol = 0.031 + 0.924] + 2.08k 
Oy = — (>.031 + 2.77) — 2.08 k 
of = 0 .9001 -r 2 . 7 7) + 4 . 07 k 
of= 12.951 + 0.924)— 0.099k 

Oh = 1 .381 4- 0.924) + 2.22k 
()f= — 1.381 + 2.77) — 2.22k 
()‘d= 5. 541 + 2.77) + 4. 20k 
( )*v — S.301 4 0.924) — 0.024 k 

();„*= 3 . 1 91 + 0 .924 ) + 0.0595 k 

()+ -3.1 91 + 2.77) - 0.0595k 

Oyh = 3.74i + 2.77)+2.04k 

OlM= 10. 1 11 + 0.924)4- 1.92k 

By using these atom position vectors, the lattice 
of A lot ) : j can be reconstructed. Projections of this 
lattice are shown in figures 8 and 9. In these 
figures the cell boundaries are shown and the ('('li- 
ters of the oxygen atoms are connected to outline 
the octahedra. Tt is readily apparent that Mo0 3 
has a layer st met ure. The layers are formed from 
one kind of octahedron and alternate octahedra 
are oriented differently. The layers are staggered 
to form a semi-interlocking joint. The three 
“kinds” of oxygen atoms result from the particular 
way that the octahedra are joined in this structure. 
The ()[ atoms are common to three octahedra and 
are a part of an edge joint. The () n atoms are 
common to two octahedra and form a corner joint. 
The On, atoms belong only to one octahedron 
and are unshared. Presumably, however, the 
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() n , atoms play the dominant role in the interlayer 
bonding. 

The bond vectors have been computed and are 
given as follows. As in the case of Mo0 2 , one 
octahedron is singled out to illustrate bonds and 
interbond angles. The bond vectors point from 
the molybdenum atom at the center of the octa- 
hedron to the oxygen atoms at the cornel's. 

In the +a direction, 

\lo iv ()^i== 1 .Sli + 0.270k (L= 1.83 A) 

in the —a direction, 

Mo^(>T=-2.27i + 0.430k (L=2.31 A) 

in the + b direction, 

Mo-OT ==— 0.4S(M+1 .S5j+0.240k (L=1.93 A) 
in the — b direction, 

MoT)T =-0.4801-1 ,X5]+0.240k (L= 1 .93 A) 
in the +c direction, 

Mc7 v OiT = 2.08k (L=2.08 A) 

and in the — c direction, 

Mo :, (fe = - 1 .891c (L— 1 .89 A) 

Vectors for the Mo— Mo bonds are directed 
from the nearest neighbor Mo atoms to the Mo iv 
atom at the center of this octahedron and are given 
as follows (the subscript X(J indicates next cell); 

In the — a + b direction, 

M o iv "\To^ — —2 . 7 (>i + 1 .85] T 0 . 070k (L=3.3S A) 
in the — a— b direction, 

Mo iv Moj“.=— 2.701— 1.85J+0 070k (L=3.38 A) 
in the -f-a + b direction, 

Mo i \\lo’ 1 -4.17i+l .85] - 1 .32k (L=4.75 A) 

in the +a— b direction, 

Mo‘^lo£ =4 . 17! — 1 .85] - 1 .32k (L=4.75 A) 

in the +c direction, 

Mo^Mo£= 3.97 k (L - 3.97 A) 


and in the — c direction, 

Afo^Mo* — —3.97k (L^ 3.97 A) 

The O — O bond vectors ar c* given for the bonds 
within an octahedron. The interlayer 0 — 0 
bonds will be considered subsequently. The se- 
lection of positive direction for these vectors is 
arbitrary and not of any significance. The O — O 
bond vectors were chosen positive toward O m as 
much as possible and positive toward 0 It other- 
wise. These vectors were obtained by deter- 
mining the vector difference between the appro- 
priate two Mo— O bond vectors as follows (the 
subscripts A, B, L, and R designate, respectively, 
atoms above, below, left, and right with respect to 
the chosen orientation): 

( )iOVT= 1 .791 — 1.85] — 0.19k (L=2.58 A) 

(>ioC— 1.791 + 1.85)— 0.19k (L=2.58 A) 

OjOu a ~ 2.2 / i + 1 .65k (Ij=2.81 A) 

f — 2 .271 — 2.32k (L=3.24 A) 

()j‘ L Oib = 2.29i + 1 .85] + 0.03 k (L=2.94 A) 

O m Oih=2.29i — 1 .85] + 0.03k (L=2.94 A) 

Oi v , A 0ib=1.81i — 1.81k (L— 2.50 A) 

OiiB0in= l.Sli + 2. 10k (L — 2.82 A) 

0|f j OiiA = 0.4Si T 1 .85] -j- 1 .84k (L— 2.05 A) 

O? R (isI=0.48i— 1.85] + 1.84k (L=2.05 A) 

OiL0nB=O.4Si + 1.85j — 2.13k (L=2.80 A) 

0¥ k ( = o .4Si - 1 .85] - 2 . 1 3 k (L = 2 .80 A) 

The interbond-angle relationships are perhaps 
more dearly demonstrated by establishing a sys- 
tem of polar angles. 'Phis system shows the nature 
of the distortion of the octahedron better than 
using ihe angles between bonds. There is a plane 
of symmetry containing the Mo lv , two Oft, Oih, and 
Oj atoms which provides for easy orientation of 
the o<* ahedron. 

A unique and important feature of the Mo0 3 
lattice is the layer structure. The O n i atoms are 
unshared by oetahedra. If the nearest neighbor 
relationships only are considered, the Om — Om 
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bonds make up the interlayer bonding. The Ojh 
atom has been arbitrarily selected as a represent- 
ative atom and the bond vectors to O n[ atoms 
in the next layers have been determined as follows: 

-0.551- 1 ,sr,j + 1 .99k (L— 2.77 A) 

b¥iiOi„, 3= -0.551-1 .S5J-1 .9Sk (I.— 2.76 A) 


OSiO]i^I^-0.r>5t + 1 .<S5j + 1 .99k (L— 2.77 A) 

6vi.<wl= —0.551+ 1 .X5j - 1 -9.sk (L = 2.76 A) 

OSiOr=-2.S4t+2.03k (L=T49 A) 

The last bond vector represents the bond between 
a layer surface atom and the nearest nonsurface 
atom in the next layer. 
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APPENDIX B 


METHODS OF ENERGY CALCULATIONS 


There are four sections in this appendix which 
are devoted to the numerical work involved in 
(he calculation of the contributions to the lattice 
energy. The theory which forms the basis for 
the practical calculations discussed herein has 
been introduced in a general way in the section 
entitled “Introductory Discussion of Theoretical 
Concepts” and the particular development of 
the equations has been given in the section 
entitled “Lattice Energy.” These equations are 
used in this appendix to calculate the lattice 
energy of molybdenum dioxide and molybdenum 
trioxide. There are four sections in which 
coulombic energy, polarization energy, Van der 
Wauls energy, and repulsion energy are discussed. 

CALCULATION OF COULOMBIC ENERGY 


The coulombic energy of a crystal is given by 
equation (24) as follows: 


l T c- 


\ av c 2 A 

2* 


In order to have the coulombic energy expressed 
in kilocalories per mole, this equation is modified 
as follows: 


where 


( * N Av e 2 A 

!()-* 2£ 


(79) 


X Av Avogadro’s number, 0.024 X 10 23 molecules/ 
mole 

e conversion factor, 2.1)889 X 10" 11 keal/erg 
molecules per cell 

e electronic charge, 4.802X10“ 10 esu 


The constants may be combined to give a working 
formula (eq. (78)): 


Pc — — 


105.9A 


kcal/mole 


The factor A is determined by equation (22) as 


follows : 


=2 2 

7J = 0 ij 


ZiZ, 

•III 


The limiting value of this equation is obtained by 
taking finite sums as follows (eq. (21))): 


k 7 7 

a,=S2 I 1 

ij = 0 ij O J j 


The (p antity A* is determined on a computer for 
values of k from 0 to 0 or 8 where the convergence 
should be very good; then, graphical means are 
used to get the limiting value A at k equal to 
infinity. Two systems of graphical extrapolation 
were used in this work. One system was to plot 
A* as ordinate and A* — A*_, or AA* as abscissa 
then extrapolate to A at AA k equal to zero. 

The other system was to select a convenient 
value if A* at a point where the sum was eon- 
verging uniformly. This starting point may be 
designated A b . The following quantities are then 
determined and plotted: 


log e (A« — A h ) is plotted on the ordinate 

1 /k is plotted on the abscissa 

The ex rapolation to 1/k equal to zero, or k equal 
infinity is made and the limiting value 
(A« — A.) limit is obtained. Then, 


A — A b T (A« — A h ) 


Both i iet hods wore used and the two extrap- 
olation? were adjusted to give values of A that 
agreed and at the same time conformed to 
reasonable curvature of the extrapolation curves. 

The most complicated part of the work in com- 
puting the coulombic energy is the computation 
of A which, in turn, requires computation of a 
scqueiue of values A*. This work must be done 
on a computer and the following comments are 
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offered to aid in the programing of the problem. 
The input data are as follows: 

(1) Position vectors of all ions in a unit cell. 
The form of these vectors is (eq. (31)) 

I j = ] jia + .I 1 ja7 

(2) The angles a, p } and 7 between the axes of 
the reference system 

(3) The cell dimensions a, b, and e in angstrom 
units 

(4) The values for l lf l 2} and / 3 determined from 
values assigned to 77 

(5) The values Zj which represent the number 
of charges on the jth ion and the sign of the charge 

Equation (23) which is 


a .=ss 

1»=0 ij 


z,z, 


represents interaction between all the ions, desig- 
nated by i, in the cell at the origin and all the 
ions, designated by j, in a cell located by R where 

(eq. (32)) 

R ~ /j a a “| _ /; _ b/3 -(- / 3 eT 


It is convenient to form a matrix 

b I- h . . . h 

T Tii 1 1 2 1 13 • • • l Ij 

T, L, I 22 1 2 , . . . l 2j 

bi T33 1 32 I33 • ■ ■ bij 


iiLiii iis 1.3 . . . i.J 

where the jtli ions are arranged across the top and 
the it h ions down the side. Each position in the 
matrix may be designated by the symbol hj. 
Then, this symbol may bo assigned to represent 
the following operation: 


_z 1 z j 
ll_ < 1 ./ 


(si) 


Terms from all positions in the matrix are added 
and this large expression is then expanded by per- 
mitting r? to increase from zero up to the point 
where the series converges. It is important to 
note that when R~() and >? = () then all terms on 
the diagonal must be discarded. Thus, 


Then, if 1 1 and Jj an* the position vectors of the 
ith ion and jth ion in their respective cells, the 
relative positions of these two ions is expressed by 
the vector (eq. (3b)) 

lhj~ I >1)37 

Tin* coefficients of the unit vectors are (eq. (37)) 
I) iji — I ji + / 1 a — 1 11 

e„ 2 — ij 2 +y>“i .2 

E j J;J — I j3 + / :i c — 1 13 

Prom this, du may be obtained as the magnitude 
of j lie vector Dij as follows: 

d u= (U jj ■ D ij) 1/J 

= (Dfji + Df J2 +D“j3^2D| jl Di J2 cos 7 

+ 2D, j 1 1) U 3 cos /3+2Dij 2 Dij3 cos a) 1/2 (SO) 


(N2) 

ij =0 ij 17=0 ij Ojj 

and the special provision for 77 — 0 must be tacitly 
understood. 

These methods were applied to two oxides of 
molybdenum. The convergence of A* is shown 
in figures 11 and 12 for Mo() 2 and Mo0 3 , 



v 

FiciUBE 1 1 . — -Convergence of A* for Mo<L. 
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V 

Fkii kk 12. (.’oiiVLTgenec of A* for M 0 O 3 . 
respectively. Tho limiting values were found to 


he: 

For Mo() 3 , 


A — 50.S, { — 4 


and for Mo() ;{ , 


A— 141 , £=4 

These values were substituted into the formula 
<eq. ( 7 S ) ) 

1 <>;>. 9 A 


to obtain eoulombie energies as follows: 
For Mo(),, 


l" r — —2,1 10 kcal/molc 
and for Mo() 3 , 

l\.= — 5,Sf)0 keal/mole 

r Th est % computations required about 30 minutes 
each on an IBM 704 electronic data processing 
machine*. 

CALCULATION OF POLARIZATION ENERCY 

The polarization energy of a. crystal is given as 
follows: 

i>--s NAv r i< ‘ i (4) ) -(A,) 1 (x3) 

where 

\ Av Avogadro’s number, 6.024X10**, molecules/ 
mole 

c (‘(inversion factor, 2.3NN9X 10“ M > keal/erg 
oj |>olarizal)ility of jth ion, cm 3 
e electronic charge, 4.S02X 10“ 1W , esu 

and the index j represents the jth polarizable* ion 
or atom in a molecule and tin* summation is over 


all ions or atoms in the molecule. The constants 
maybe* *ombined to give* a working formula, in kcal / 
mole, 

0 s4 ) 

j 

where tho prime indicates that the polarizability 

is in A 3 units. The (‘valuation of the vector (a p ^ 
must be made 011 a computer. The following 
comments are offered to aid in the programing of 
the problem. The input data are as follows: 

(1) Position vectoi’s of all ions in a unit cell; 
the form of these vectors is (eq. (31)) 

Ii^I n a + liJ+Ii37 

(2) Position vectors of the selected representa- 
tive*. polarizable ions. These* will he designated 

G>. One ion of each kind j will be* selected 
— ^ 

from ions I,. All the* ions j will constitute a 
molecular unit. The* form of these vectors is 


fij.iT (N5) 

(.'?) Tl)(> ilMglt'S a, fi, :lll(l 7 l)(‘( WIR'D till' it.Xl'S of 
t lie* red ‘renee* system 

(4) Tin* <*(*11 dime'iisioiis a, b, and e* in angstrom 
units 

(5) The values for l> , and / 3 determined from 

values issigm*d to 77 (in t his work only values 0 and 
1 are assigned to 17) 

(0) 'The value's Zj which represent the number of 
charge; on the it h ion and tlu* sign of tin* charge* 

Tin* vector (a,,), is computed from equation 
(40) w lich is 



where cep (MO)) 

l).j=I>i.ia ■ D.Jk I>i J3 7 

and (eqs. (37)) 

D iji 1 1 1 11 I pi 

— gj> — 1 12 — ^->h 
0 1^— gj 3 I ta /3O 
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The expanded expression for 




= W T — - 

«j=n i (Dfj! +D£2+D- ! j3 +2D 1 jiDij2 cos 7 


Z|(DijiO! + D l}2 /3 “b D{j 3 t) 


+2D| ja Di J 3COsa+2D t j 1 Di J3 cos p) 3r2 


(Sh) 


In this expression lor 77 — 0, there will he terms 
—> - » 

where Gj = ] , and these must be discarded to avoid 
dividing bv zero. They would represent the 
interaction of an ion with itself. 

In computing' the polarization energies of the 
two oxides of molybdenum, the polarizable ions 
were select ed as follows: 

For Mo() 2 , 

(V ^ lo^ - 1 .295a -)- 2.42 1 3 + 2.899-y 

G 2 = (>i = Odil 4a + 1.017/8 + FTitvy 

^3 = ( >n = 2. 1 78a -j- 3.389(5-!- 1 .(>827 
and for Mo() 3 , 

G, = ^ 1 o*" ^ 5 . 54 a + 2 . 7 7 j8 + 2 . 32 7 
G, = ( )i = 0 .03a + 0.924/8 + 2 .O87 


and 


- - 0.680a — 0 .00686/9 — 0.07 1 \)y 


and for M0O3, 

(aJ.^-o.cs 

( A p ) 2 = -0.:mS —0.0625/3—0.2827 
(a„) 3 = — 0.853a — 0.0031 1/8+0.2897 
(4)4 = — 1 ,295a — 0.00362/3— 0.09057 

(a p ),-(a p )i = 0.4676 (Mo 6 +) 


(-0,(4) 


2=0.1943 


(O*-) 


(4)3.(4)3=0.8111 

(4)4 (4)4=1.298 J 

Tlie polarizahilit ies used w(*re as follows: 


G, ^ Oj£ 5.54a + 2.77/3 + 4 . 2 O 7 
Ch—Oih — 3.74 a -*2.77/3 -*-2. 047 

Tlit' coinput oil CO, vectors were ns follows: 
For Mo() 2 , 

(a„), = 0.304a-0.0473S+ 0.05997 
C0., = 0.242a 4 0. 1 73/8+ 0.3537 
(a p )s“ 0.07 13a — 0.0342/8 — 0.002627 

and 

( A , , ) 1 - ( A p ) , = ( > . 0 7 0 r> (Mo‘+) 

(-4) 2 -(4X- = 0.12r,3 ) 

(O 2 -) 

(a i> ) 3 -(a I) ) 3 == 0.00<i4r>) 


For Mo 4 +, 
for Mo" + , 
mul for O 2 " 


a' = 0.35 A 3 
a' = 0.26 A 3 
a' = 2.22 A 3 


These (plant it ies may he subst it uted into t he work- 
ing formula (eq. (84)) 

U p = - 16o.92&'(4)r(4)j 

j 

to yield polarizat ion energy. Tin* results are 
for MoOo, 

Up — 53. 1 2 kcal/inole 
and for Mo() 3 , 

l ,, = 808.40 keal/mole 

Aj)j)roximately 10 minutes on tin* IBM 704 com- 
puter were required to compute a set of electric- 
field vectors for either oxide. 
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CALCULATION OF VAN HER WAALS ENERGY 

The Van tier W mils forces are very short -order 
forces and consequently, the computation of the 
Van der Waals energies may be limited to con- 
sideration of nearest neighbor interactions only. 
In computing the Van der Waals energies on a per- 
mole basis, it is necessary to sum the contributions 
of all bonds associated with one molecular unit 
and multiply t his by Avogadro’s number. F or one 
bond, the dipole-dipole and dipole-quadrupole 
energies are given by the following expressions: 

For dipole-dipole energy (eq. (62)), 


must >e estimated since it does not have a closed- 
shell structure. It is assumed that the two 
d-elee runs will cause an increase in X for Mo 4 ' 1 ' 
over that for Mo fii ' and a value X = 5.5 is estimat- 
ed for Mo 4+ . 

In order to establish a working formula for the 
dipole -dipole energy, equation (54) must be mul- 
tiplied by e 2 a 0 I/2 which converts from atomic to 
cgs units and where 

e electronic charge, 4.802 X 10 -10 esu 
a 0 Bohr radius, 0.529X10“* cm 

The following constants must also he applied in 
order to give a final expression in the units 
keal/i Hole: 


and for dipole-quadrupole energy (eq. (04)), 

/r T \ 

tVdtiMjk 8 


In these expressions the indices j and k refer to in- 
teraction between an ion of element j and an ion of 
element k. The index i refers to bond distance i 
for those two (dements. In the distorted struc- 
tures there are several bond distances for each pair 
of (dements rather than a single distance found in 
more symmetrical structures. The coefficients 
J lk and L jk are determined as follows. The dipole- 
dipole coefficient is computed by using equation 
(54) which is 


where 


" Jk~ 


a j<*k 


/*J . 

V.N, + 


/ « k 

VXk 


a polarizability 

X number of polarizable electrons in the ion 
For ions with closed-shell structures, the value 
for X T is that computed for the isoeleetronic inert 
gas. These values were computed by TIerzfeld 
and Wolf (ref. 13) as follows: 

For argon, 

X = 4.6 

for krypton, 

N — 5.1 (inter])ola ted) 
and for xenon, 

A — 5.6 


On this basis, the value of X for Mo 6+ is 5.1 and 
for O 2- is 4.6. The value of A for the Mo 4+ ion 


c 2.3889 X 10“ 11 kca 1/erg 

X Av Avogudro’s number, 6.024X16“* mole- 
cules/inole 

These constants arc* combined to give th(* formula, 
in kciil/mole. 


3 


(I (Id)ljk - ~24 1 .4 1 / y 




/< 


\ W } Vn 


«k 


(«) 


it is convenient to use the form, in kcal mole, 


(L T dd)ijk— — 24 1 .4 1 


when 


C*= 


la’ lal 

V-n+Vn. 


Th * polarizabilities a are in A 3 units and Ri is 
in A units. The dipole-quadrupole coefficient is 
computed by using equation (57) which is 


1 


— '} ( Jkk \ 

Jjk c 2 \A jtq X\a k / 


The \ nit conversion factors are combined to give 
the following formula, in kcal/mole, 


when 


(Udu)uk- -4,060 R Jk 

,i _ J! f i. w, 

Jk 4.80* \N .a. Stfitv) 


(90) 


( 91 ) 


*\Nja, N k a k v 
The final formulas for computing the Van der 
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Waals energies, in keal/mole, are as follows: 


U«=-241.41S :M ™? J ' 

ijk n. ] 


u rt<1 =- 


-4,060 X! 

ijk 


H." 


(92) 


(«w) 


where M 1Jk is the number of bonds present between 
ions j and k with a distance Hi. The terms in 
the summations of equations (92) and (93) are 
tabulated in tables 1 and 2. These values are 
substituted into equations (92) and (93) to yield 
the following Van dor Waals energies, in keal/mole: 


For M 0 O 2 , 


and for Mo() 3 


l T dd =- 105.99 
l T d<j— — 7(>.57 


Udd = 

IV 


— 72.70 
-25.13 


CALCULATION OF REPULSION ENERGY 

The repulsion energy l r is given by equation 
(4) which is 

Ur — — Pr(l e + 41-p+bl rtrt + Sl ^,) 

where the U terms are the different kinds of 
attractive energy and p r is a repulsion coefficient. 
The repulsion coefficient is given by equation 
{74 ) which is 

(i’i + 1 

<iu n,j 


SM„Ar 

ij 


XI MijUijftj 

ij 




where («*<(. (66)) 


Aj = l + |Y+ 


Z, 

I’j 


and 

Mij number of bonds per mole between ions i 
and j at a distance d tj 

r ionic radius 

dij bond distance between ions i and j 

Z valence of the ion 

P number of outer electrons for the ion 
Hi j a constant dependent on the particular ion 
pair 

The value of n tJ is obtained by using the method 
of Pauling (ref. 2, p. 509). The quantities used 
in the calculations of p r are given in tables 3 and 
4. The calculated values of p r are as follows: 

For M 0 O 2 , 

p r — 0.1 295 

and for Mo() 3 , 

p r = 0.12(H) 

Substituting these values in equation (20) gives 
the repulsion energies ns follows: 

For M 0 O 2 , 


U r =4(i2.44 keal/mole 
and for Mo() 3t 


l T r = 1,201.12 keal/mole 



APPENDIX C 

METHOD OF OBTAINING POLARIZABILITIES OF IONS 


Tin* polarizabilities of the ions assumed to exist 
in the oxides of molybdenum must be known in 
order to compute lattice energy. Since some of 
these data were not found in the literature it was 
necessary to search for an empirical relationship 
that would be capable of providing reasonably 
accurate values of the polarizabilities. A method 
frequently employed is the use of l lie cube of the 
ionic radius as the polarizability. This value is 
generally recognized as a rough estimate; however, 
it does not agree with the method established in 
this investigation. The empirical relationship 
presented here is believed to be a new and more 
accurate method of obtaining polarizabilities of 
positively charged ions. 

The polarizabilities of () 2 ~, Mo 41 ', and Mo ,i + 
ions are needed to compute the polarization energy 
and the Van dor Wauls energy. The polariza- 
bility of an ion is sometimes regarded as a con- 
stant. Actually, tin 1 polarizability of an ion varies 
over a small range of values according to the 
environment. Ionic polarizabilities have been 
calculated and measured by many investigators. 
A small part of this work is described in references 
14 to lb. 

According to reference 15 the polarizability of 
varies according to crystal type. Of the 
available data, it appears that Ti() 2 provides an 
ion environment most like that in the oxides of 
molybdenum for the () 2 ~ ion. Measurements on 
TiO* yield a value of 2.22 A 3 for the polarizability 
of the O' ion. Fa j a ns and Joos (ref. lb) found a 
value of 2.75 A 3 for the () 2_ ion in solution. 
Polarizability calculations for the () 2 ~ ion in a free 
or gaseous environment yield values near 5.8 A 3 . 
Polarizability measurements on several crystals 
yield values ranging from 2.0 to 2.5 A 3 for the (F - 
ion. The value of 2.22 A 3 as tin* polarizability 
of the 0‘ - ion has been adopted for the computa- 
tions in this paper. 

The polarizabilities of Mo 4+ and Mo fi+ ions were 
not found to be reported in the literature and 


presumably have not been determined. An em- 
pirical relationship was determined to predict the 
polarizability of ions having the inert-gas-type 
closed-outer-shell structure. It is found that the 
following relationship holds for ions having the 
same charge and closed-shell structure: 

a'-k.v (84) 

where 

a polarizability in A 3 

r ionic radius in A 

Ej ionization potential in electron volts 

k a constant for ions of the same charge 

Polarizabilities of ions in solution measured by 
Fajans and Joos (ref. lb) are plotted against 
r'VE! in figures 15 to lb with a separate plot for 
ions of different charge. The slope of each 
straight line produced yields a value of k for the 
partie liar ionic charge. Values of k corresponding 
to charges of +1 to +4 were determined. A plot 
of k against charge (fig. 17) was then made and k 
was extrapolated to give a value for ions of 
charge 4-b. By this means it was possible to 
arrive at a value for the polarizability of the Mo ftt " 
ion. Hie values of k were determined as follows: 


( Miar^c | k 


y i 

12.5 

- 2 

24 

- 5 

44 

■ 1" 4 ; 

70 

1 \ 0 1 

tool „ 

j-ti 

1 44 j * <lX trapol:it ion 


For Mo* + , vv =0.00192 which yields polarizability 
a equal to 0.2b. 


to 
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Ion charge 

Fici kk 17. — Polarizability constant for cations. 

Presumably, the polarizability for the Mo ,f ion 
ran be obtained only by estimation. There are 
some factors to serve as a qualitative guide in such 
estimation, however. The polarizability of the 
Mo 4f ion must be larger than that of the Mo wi ion. 


The Mo fi+ ion has eight electrons in an outer- 
closed-shell configuration which is inherently 
stable. The Mo 4+ ion has the eight electrons in 
the closed shell plus two additional electrons. 
The additional two electrons are in larger orbitals 
and are less firmly bound. They would tend to 
cause increased polarizability. However, the 
presence of the two outer electrons would inhibit 
somewhat the deformabilily of the eight electron 
orbitals in the next inner shell and reduce the 
polarizability slightly. In consideration of these 
relationships, the polarizability of t ho Mo 4f ion is 
estimated to be 0.115 Ah Xo doubt there is a 
good possibility of error in arriving at a value of 
the polarizability of Mo 4+ ion by this means. 
There was no other apparent recourse and it is 
felt tint the error introduced is less than that 
which would result from assigning a zero value to 
the polarizability of the Mo 4+ ion, that is, from 
ignoring it altogether. 

In summary, the polarizabilities used in com- 
pulations in this paper are as follows: 

For 0~x 

a' = 2.22 A 3 

for Mo< 1} , 

a' = 0.20 A 3 

and for Mo 4+ , 

a' = 0.:i5 A 3 



APPENDIX D 


THERMOCHEMICAL DATA AND THE BORN-HABER CYCLE 


The values which have been obtained for the 
lattice energies of the oxides of molybdenum by 
theoretical methods may be compared with experi- 
mental measurements. Thermochemical meas- 
urements are made of the heats of formation of 
compounds. Idlest* data are transformed into 
lattice energy by the Born- Haber cycle. The 
thermochemical data are given in the first section 
of tins appendix and these are applied in the 
second section which describes the Born-1 label 1 
cycle calculations. 

THERMOCHEMICAL DATA 

The heats of formation of Mo(b and MoO d have 
been measured (refs. 17 and IS). Tin* recom- 
mended mean values are given as follows (at 
298.15° C): 

Mo(c)t () 2 (g) >Mo() 2 (c) AH— — 140. S 

kcal/mole 

Mo(c)+ lUOdg) >Mo(),(e) All= — 17S.1 

k(*al /mole 

Other thermochemical data from reference 19 are 
Mo(e) >Mo(g) AH — 155.5 

{X, >20 All — 1 10.8 

Electron Affinity: 

Tin* (*lectron affinity of () 2 “ has been recently 
determined by Huggins and Sakamoto (r(*f. 8, p. 
241) as follows: 

Electron affinity O 2 " 102±15 kcal/mole 

lonizat ion Potent ial: 

The energies required to remove the first six 
electrons from molybdenum are given as follows 
(ref. 20): 


For 

Mo, 





I_ . 

.7.18 

electron volts 

IV 

40 

.5 electron 






volts 

ir 

15.2 

electron volts 

V . 

50 

electron volts 

m 

27.0 

electron volts 

VI 

72 

elect ron volts 

for 

Mo 4+ , 






89.88 electron volts or 2,(108 kcal/mole 
and for \Io 6+ , 


217.88 electron volts or 5,014 kcal/mole 

LATTICE-KNKROY CALCULATIONS BY BORN-HABER CYCLE 


The Born-Haber cycle is (*x|)ressed by the 
following equat ion: 



AH r - AH f -AH K -AH d — E a-E, 

(95) 

wl 1 (']'(> 
AH, 

lattice enthalpy (^lattice energy) 


AH k 

fieat of atomization (molybdenum) 


All, 

heat of dissociation (oxygen) 


K, 

electron affinity (oxygen) 


K, 

ionization potential (molybdenum) 


AHf 

heat of formation 


Also, 

< 

<1 

li 

< 

(90) 


The term AnRT is relatively small and is 
neglected; therefore, AH l «AE l . 

For MoOo, 

AH l — - 140.8- 1 55.5- 1 1 0.8 — 524 — 2,008 

— 2,805 kcal/mole 

For Mo() ;i , 

Allj,— — 178.1 — 1 55.5 -175.2 — 480 — 5 ,0 1 4 

— —0,009 kcal/mole 
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